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SYNOPSIS 

Vinod Tare 
Ph.D, 

Indian Institute of Technology, Kanpur 
September 1981 

MODELLING AND SIMULATION OF DEEP- BED 
FILTERS FOR '.^ATER CLARIFICATION 

Filtration is a subject which has, in recent years, developed 
a rich and varied li-terature. Despite the dramatic progress in 
understanding filtration process, a canaprehensive method capable 
of predicting the dynamic behaviour of a filter over its en-tire 
operational span has remained out of reach. No doubt, ■t±ie 
uncertainties of the nature of the deposition process, particu- 
larly the critical ones relating to the nature of -the distrilxation 
of -the deposited matter wi-thin -the media and its evolution wi-th 
progressive deposition, have stood as the main impediments obstru- 
cting progress towards such a general predictive scheme. In the 
present research this problem is resolved by postulating an over- 
all picture of the filtration process which incorporates -the 
significant role pUayed by the retained particles on filter 
performance as the filtration progresses, 

Ihe ultimate objective of the theoretical approach is to 
derive mathematical formulations which quantitatively describe 
■the conpllete c^namics of deep bed filtration, taking into accxunt 
the compHexed interaction be-tween the filter grains and liquid 
solution partied. es to be removed as well as the effects of 
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operating and syston variables, At this stage of its development, 
it has been possible to predict with reasonable accuracy the head 
loss and efficiency of cxLean filters. During filtration, both 
head loss and removal efficiency change with time. These changes 
are caused by particles that are retained in the filter bed. In 
the present research an attempt has beon made to formulate a new 
conceptual model for predicting the filter peirformance as the 
filtration progresses based on the concept that the filtration of 
suspended particles in deep bed filters is a two stage process, 

Ihe first stage involves the transport of suspended particles to 
the media grain or already deposited partides and the second 
stage acccunts for the attachment of these particles to the media 
grain or the d^osited particles. It is postulated that the 
deposited partides act as collectors, and some fraction of the 
deposited partides (Contribute for the reduction in effective 
area of the media grain and retained partides acting as collectors. 
To provide a physically realistic model as basis for quanti- 
tative descriptions of the* various processes of interest which 
occur in the xcorous medium, the concept of Unit Bed Element (UBE) 
and of the unit collector is used. The filter bed is postulated 
to consist of a number of unit bed elements in series. Each one 
in turn is assumed to consist of a number of single collectors. 

The capacity of the bed to remove particulate matter can be 
described by the collection efficiency of the unit bed element, 
(lefined as the fracticxi of susp^ided particles entering the UBE 


removed 



At the first stage of model developrienty the filtration 
efficiency is described by a set of three interd^endent equations 
for monosized particles in suspeision assuming the porosity of the 
bed to ranain constant as the filtration progresses. The first 
equation computes the single collector efficiency related to clean 
grain efficiency and retained particles acting as collectors. The 
second equation estimates the number of retained particles acting 
as collectors as a function of depth, time and local particle 
concentration in suspension, The third equation predicts suspended 
particle concentration with time over the entire filter depth. The 
head loss development is predicted ty modifying the classical head 
loss equation for clean filters incorporating change in the 
surface area to volume ratio of the filter grains. 

In the second stage, the model is refined to eliminate and 
check the assumption of constant porosity . The change in porosity 
is computed from the decreased pore volume because of the retained 
particles. Finally, the model is extended to include the effect 
of multisized particles in suspension as well as the multisized 
media. 

The model is verified with a set of laboratoay data obtained 
from literature and the experimental investigation carried out • 
during the present research. Uniform latex particles with mean 
diameters of 0,109, 1,1 and 7,6 microns were used as the turbidity 
in the fil-ter influent in most of the experiments. Appropriate 
experiments were conducted in which the size and concentration of 
suspended particles were varied while physical and hydraulic 
filtratibn paraneters were maintained at constant levels. The 
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model predicted values very well agree with almost all experimental 
observations. 

The simulation of deep bed filters is achieved in two steps. 
The first step involves the simulation of filter performance 
characterized by the effluent quality and head loss development. 

The second step includes the simulation of filter state character- 
ized the clogging zone. The clogging zone is obtained by 
proposing an equation for the clogging front. The filter perfor- 
mance and filter state are simulated for various physical and 
physico-chonical parameters influencing water and wastewater 
filtratison such as size, density and concentration of suspended 
particles, size of the media, initial bed porosity, flow rate, 
particle- to-media grain attachment, and particle-to-particle 
attachment by experimenting on the proposed model. The model is 
able to predict the filter performance in accordance with the 
observations made in laboratory and field. 



1 . INmODUCTION 


Filtraticsn, a physico-chemical process is one of the 
most important unit operations anployed for the final clari- 
fication of water and wastewater in treatment plants. The 
process involves the retention of suspended and colloidal 
particles by the filter grains and/or previously retained 
particles as the suspension is made to pass through the 
granular filter bed. The main indicators of the dynamic 
behavior of the operation are filtrate quality history and 
the pressure drop history requir«3 to maintain a uniform 
throughput. The first defines the effectiveness and the 
capacity of the filter bed, and the second determines the 
length of the filter run. 

The importance of an accurate simulation of the dynamics 
of packed bed filtration is evident, as it provides the 
basis for rational design and optimization. Recent advances 
in the study of deep bed filtration have shown that inspite 
of the inherent complexities of the process, the problem can 
be studied fruitfully based on fundamental consideration. 
Resurgent research efforts of the past two decades in this 
relatively old engineering practice have already yielded 
important results of practical significance. In general, 
the study of deep bed filtration has been made using two 
different yet often ccanpiementary approaches: phenomenological 
and theoretical. The phenomenological approach describes 



the dynamic behavior o£ deep bed filters with the use of a 
set of partial differential equations and characterizes the 
filtration mechanisms by means of several model parameters 
the values of which ^ for a given application, can be obtained 
from appropriate bench scale experiments. Ihe solution of 
these equations provides the basis of design, scale up and 
optimization. Although this approach affords only limited 
insight about the physical process of filtration and is 
therefore not entirely satisfactory from scientific viewpoints, 
it represents a practical and rational methodology to the 
design of industrial scale depth filters, provided the 
bench scale experiments frm which the rel event parameters 
are obtained can be conducted with sufficient accuracy. 

Ihe ultimate objective of the second approach viz,, 
theoretical, is to derive mathematical formulations which 
quantitatively describe the con^lete dynamics of deep bed 
filtration, taking into account the complexed interaction 
between the filter grains and liquid solution particles to 
be removed as well as the effects of operating and system 
variables. At this stage it has been possible to predict 
with reasonable accuracy the performance of clean bed filters 
in terms of head loss developnent and removal efficiency . 

The loss of head at the start of filtration can be predicted 
by knowing the media size, filtration rate, and bed porosity. 
Similarly, the removal efficiency of clean bed filters 
depends on the characteristics of suspended particles and 
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filter media like size, density and surface properties, the 
filtration velocity and the depth and porosity of the filter 
bed , 

During filtration, both head loss and removal efficiency 
change v/ith time. These changes are caused by particles that 
are retained in the filter bed. O'Melia and All (1978) 
proposed a model for the developnent of head loss and removal 
efficiency during filtration by packed beds based on the 
postulate that some retained particles can act as collectors 
and thereby improve filtration efficiency. Though, the 
concept seems to be appropriate and the mathematical model 
agrees within the certain limits of experimental data, the 
estimated coefficients may not reflect or may not have any 
correlation with the actual properties of t^he suspension and 
filter media responsible for filtration. This is because 
the aforementioned model fails to properly incorporate few 
important factors, such as surface coverage of the filter 
media grains and that of retained particles acting as 
collectors, variation of retained particles acting as 
collectors along the filter depth etc., which would influence 
the filter performance to a large extent. 

Therefore the present research was undertaken with the 
main objective to formulate a new conceptual model incorpo- 
rating above-mentioned factors which would help in dynamic 
simulation of filter performance. The model should be able 
to predict the filter performance in accordance with the 
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observations made in laboratory and field. Using the 
proposed model an erjiation is developed for the depth of 
floe penetration which vrould help in evaluating the optimum 
filter depth. 
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2. LITERATURE REVIEW 

Filtration of water through a porous media is an old 
process in the field of water purification. Slow sand filters 
were introduced in Scotland in 1804 and into the United States 
in 1872 (Fuller, 1933), Rapid sand filters came into existence 
towards the end of the nineteenth century and since then it 
has become a major unit operation in water purification plants. 
The main function of rapid sand filtration is the final clari- 
fication of settled water. Many attonpts have been made to 
enhance this achievement in order to cope with demands for 
larger amounts of higher quality water. Thus a considerable 
interest in the practice and theory of water filtration has 
been maintained. Early developanents were focused on refining 
the structural features of filter; this was followed by a 
period of refinement in filter operation. Subsequently, 
emphasis was placed on proper conditioning of the water prior 
to the filtration. In recent decades attention has been 
directed towards the understanding of the mechanisms of the 
filtration process so that additional improvements and 
innovations can be obtained. 

This continuous interest on the subject of water f liter- 
ati on has generated a rather extensive literature. This 
literature is so broad that a ccanprehensive review of all 
these publications would be quite lengthy and beyond the 
scope of this dissertation. Furthermore, such reviews are 
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available elsewhere (Agrawal, 1966; Habibian, 1971; Ali, 1977; 
Tien and Payatakes, 1979), However^ to facilitate the discu-* 
ssion of the scope and results of this thesis ^ a summary of 
pertinent literature is presented here. 

In general, the study of deep bed filtration has been 
made using two different yet often complanentary approaches; 
phenomenological and theoretical. The phenomenological 
approach views the filtration process at macroscopic level and 
is aimed at ( 1) predicting the dynamic behavior of filtration 
process; ( 2) development of methodology and technique for 
design, calculation, and optimization. The theoretical 
approach studies the deep bed filtration at microscopic level 
and is intended to provide information and insight about the 
mechanisms of particle deposition, the conditions under which 
this occurs, and the effect of deposition on the structure 
and characteristics of filter grains, 

2.1 Phenomenological Approach 

2.1.1 Prediction of Filtrate Quality 

The fundamental equations describing the particle 
retention of a filter bed are the macroscopic conservation 
equation and the rate equation. The conservation equation 
is similar to that used in any fixed bed process involving 
the transport of matter from the mobile to the stationary 
phase. For an axial flow filter with a constant cross- 
sectional area, this is written as 
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Here, u is the superficial velocity, c is the volume concen- 
tration of particulate matter in liquid, z is the axial distance, 
t is the time, a is the specific deposit, and F is the poro- 
sity of filter bed. 

The assumptions involved are one dimensional plug flow 
and negligible axial dispersion of particulates. The porosity 
of the bed F changes with time as particle accumulation within 
the bed increases. If the deposited matter forms relatively 
smooth coating on the external surface of the filter grains, 

F and a can be related by the simple expression 



( 2 . 2 ) 


Here, F^ is the porosity of clean bed and F^ is the porosity 
of deposits. 

On the other hand, if the deposit morphology is largely 
of the blocking type, a different expression relating F, a, 
and the interstitial velocity results. 

Similar to all fixed bed processes, it is more meaningful 
to employ a corrected time variable © defined as 

e = ' t - f — dz (2.3) 

_ u 


In terms of the variables (z, @), Eq. (2,1) becomes 


U(l£) . (-^^) 


= o 


(2.4) 
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The difference between t and © is usually small. However, 
it may become important in the interpretation of data from 
small experimental filters. 

To complete the description, the rate expression of 
filtration needs to be specified. In general, this can be 
expressed as 




U(|f) 

3 z 


© 


= G ( Y/ c, e) 


(2.5) 


Here, t is the parameter vector of function G(y , c, ff ) . 

Based on slow sand filter data, Iwasaki (1937) proposed 
the expression 


(— ) 

^ az^ 


X c 


(2.6a) 


Here, ^ is the filter coefficient. 

If one ignores the difference between t and ©, combining 
Eqs, (2.5) and (2.6a) yields 


3 a 


= u Xc 


( 2.6b) 


In other words, use of Iwasaki* s assumption is tantamount to 
the use of a f irst-C/rder expression for the rate of filtra=. 
tion. Later, Ives ( 1960) applied Eq. (2,6a) to rapid sand 
filters successfully. The validity of this expression was 
also shown experimentally by Ison and Ives ( 19 69). 

From physical considerations one would expect ^ to 
be a function of the state of the media. This is confirmed 
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by expejriinental observations. To account for the variation 
of ^ due to particle accumulation in the bed, one may write 


X = x^(x) f^ (aS o) (2.7) 

with 

f ( a ' , O) = 1 


Here, is the value of x when the bed is clean, X is the 
parameter vector ( including quantities such as interstitial 
velocity, grain size, density difference, etc.) on which X 

o 

depends, and f^ C a', a) is the function describing the effect 
of particle deposition on filter coefficient. 

The empirical determination of the functional form of 
( a' , a) from experimental data has been attempted by a 
number of investigators. Extensive reviews on this can be 
found in Agrawal ( 1966), Ives ( 1969, 1971), and Harzig et.al. 
( 1970), Table 2,1 lists some of the expressions proposed. 
Generally speaking, most of the proposed expressions can be 
derived from a general equation proposed by Ives ( 1967, 1969) 
and Mohanka ( 1969) . 


( 1 + 



( 1 - |-) ^ ( 1 
o 


^max 


( 2 . 8 ) 


Here, a., a-, ate onpirical constants, and a ..is the 
maximum value of a to bo achieved in filtration. 

It is implicitly assumed that filter media will ultimately 
become nonretontive owing to the change 'of surface character- 
istics, increased roentrainment, etc. 
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Taole 2.1. Some analytical' solutions of filtration equations 
(Adopted from Tien and Payatakes, 1979) 


S. 

No 


Investigator 


Rate expression 


Solution for filtrate 
quality 


Iwasaki 
( 1937) 


5 c 
Sz 


X c 
o 


X Z 


= c 


Mints (1951) ||= Xc^ga 


- X z , . ,.n-l 

c. -T ® (n - 1)1 ^ 

u n=l 


= I 


T e 
n 


•ae 


T = T 


, ^ \ n=. 2 

,( a ©)_, 


n n-1 (n - 2)i 


= e 


a© 


Shekhtman 
( 1961) 


3c 

3z 


X^d ~ 


j Gr )c 


X lIuc © 
o-* O 


c \ z \ due © 

o ^o o 

G + e - 1 


4 Ives ( 19 63) 


l£ 

3z 


=. (a ^ 
b a^) c 


C 

c 


Be 


=*az 


(B^ 1 + 


2uc©(ab)^^^ 

R _ e ... . . + A 

2uc©(ab)^'^^ 
e ° ^ 
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A different approach was adopted by Mints (1951) and 
later by Litwiniszyn ( 1967) . The rate expression is assumed 
to be 


( 1 £) 


X c 


ag 

u 


(2.9) 


Here/ a is the coefficient of filtration rate expression. 

The first term of the equation represents particle deposition 
as before/ whereas the second term is introduced to account 
for the scouring of deposits/ evidence of which is given by 
the observation that in a shallow filter/ the effluent was 
found to contain particle aggregates significantly greater 
than those present in the influent stream (Mints et.^., 1967). 
Mints also argued that improved filtration efficiency/ due to 
the application of polyelectrolytes / can bo attributed in part 
to the fact that polyelectrolytes increase the bonding strength 
among the dejposits/ therefore reducing the extent of scouring. 

The shortcoming of these expressions is that all the 
model parameters can be determined only empirically. Some 
attempts have been made to develop generalized correlations 
for Xq* For example/ Ives ( 1971) suggested the following 


form of correlation; 

d. e . e 2/p - Va^g e pu d e 

. rdg^ ^eit ^iapUdg^ ^ 9n u ^ ^ ^ 

( 2 . 10 ) 

Here/ d is the grain diameter/ P- is a coefficient/ d^ is 
g 1 p 

the suspended particle diameter/ k is the Boltzmann’s constant/ 
T is the absolute temperature/ |j, is the viscosity/ 
e^, ®4 exponents, a^ is the 



12 


particle radius^ is the density of suspended particle/ ^ is 
the density of liquid, and g is the acceleration due to gravity. 
However, there do not exist sufficient data or, more importantly 
accurate data to establish such a correlation. Rajagopalan 
and Tien ( 1977) also pointed out tlxe lack of theoretical 
validity for an expression such as that of Eq, (2.10). Accor= 
dingly, the prediction of filter coefficients, as well as of 
the rel event functional expressions, must largely rely upon 
theoretical means. 


2.1,2 Prediction of Pressure Drop 


For a clean filter bed, the pressure drop versus 
flow rate relationship is goveimed by the Carman-Kozeny equation, 
that is 



150 




( 2 . 11 ) 


To account for the change in pressure gradient due to filter 
clogging, either one of the following expressions may be used; 


or 


(3 p/9z) 

( 9p/ ^z ) ^ 


= P , a) 


(2.12a) 


( 3p/3z) 

i 3p/3z)^ 


1 + F2(P . cr) 


Here, p is the parameter vector of function F^( p , c) and 
F 2 ( P / <^) / p is tte pressure. 


(2.12b) 
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A number of specific functional forms have been proposed 
for ^2 . ^ Pienvichitr, 1965; Ives, 1969; Herzig 

^ 1970), Ives and Pienvichitr ( 1965) pointed out that 
most of these expressions can be derived from a general form: 


F,( P, c) = 


( 1 — ^2 


— 1 + ^ 2 ^"^ 


2-1 + (P^ + 1) 


- q2 2 

^ t ^ • • • * 


( 2.13) 


Here, and components of p 


Although theoretical expressions similar to that of Eq, (2.13) 
have been obtained, these expressions are based on specific 
deposit morphology whidi may or may not be valid for the problem 
at hand. In practice the coefficients (p^^ and example) 

have to be determined from experimental data. 

If the specific functional form of is known, and 
furthermore if a^z, ©) is known from the solution of the conser- 
vation equation [that is, Eq, ( 2,4)] and the rate expression 
[Eq, (2,5)] , the pressure distribution can be obtained fix)m a 
simple integration as follows; 




o 



[P/ a ( z* 


, ©) ] dz‘ 


( 2.14) 


2,2 Theoretical Approach 

The phenomenological equations discussed above describe 
the dynamic behavior of deep bed filtration in terms of model 
parameters and functions. The importance of estimating these 
quantities with sufficient accuracy is rather obvious. Since 
no generalized correlations for the estimation of these 
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quantities have been established, one has to resort to tiieore- 
tical calculation. Furthermore, theoretical analysis has the 
added advantage of providing a more fundamental understanding 
of the deposition phenomenon. 

The compliexity of the filtration process has led some to 
believe that a complete theoretical fonrulation of filtration 
may not be attainable at all, Mintz ( 1966) , reporting at the 
1966 International Water Supply Congress, Barcelona, stated 
that: 


"There are complicated interrelationships among many 
factors affecting the performance of rapid filters. These 
interrelationships may vary with seasonal changes in the 
quality of raw water, chemical treatment, output, and load 
changes. Therefore, it is apparent that an attempt to work 
out an exact mathematical description, with theoretical 
constants, of the filtraticn process to hold for any condi- 
tions of filter operation is bound to fail," 

Despite of this generally accepted pessimistic view, it 
has been suggested by O'Melia ( 1965) that a reasonable under- 
standing of the basic mechanisms involved in filtration process 
may be obtained by envisioning this process as involving two 
distinct steps: (1) the transport of the suspaided particles 
to the immediate vicinity of the fil er grains, and ( 2) the 
attachment of particles to the filter grains or to the particles 
already deposited. This saitie approach has been successfully 
eirployed in the formulation of the coagulation process. 

The transport step is a physical process and is controlled 
by physical filtration parameters such as flow rate, grain size 
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and shape, bed porosity and depth, and the shape and density 
o£ the suspended particles. Particle attachment, however, is 
a colloid- chemical process and is affected by a variety of both 
ehaaical and physical parameters. Ihus by considering both 
transport and attachment steps the effects of both physical 
and chemical variables can be incorporated in developing a 
combined physical “chemical filtration model. 

2,2,1 Ihe Transport Step 

The principle involved in the theoretical calcn- 
lations in quantifying the extent of transport of suspended 
particles in a filter bed is as follows, A filter bed can be 
viewed to be an assembly of particle collectors. Using the 
terminology of Payatakes et,^. ( 1973) . a filter bed can be 
considered to bo a series of Unit Bed Elements (UBE), each of 
which, in turn, is composed of a number of collectors. Since 
the thickness of a unit bed elenaent AL is alv;ays small, Eq, 

( 2,6a) can be applied to a UBE with the assumption that the 
filter coefficient remains constant throughout the element. 
Considering the i^ UBE, let c. - and c. be the concentration 
of the influent and effluent streams, respectively, the following 
result is obtained: 

c. , 

In ■ = X A L ( 2.15) 

Similarly, the capacity of a UBE to remove particulates 
from the fluid stream is characterized by its collection 
efficiency Tj* / defined by 
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( 2.16) 


( 2.17) 


Specification of geometrical configurations for the 
collectors present in a UBE is admittedly somewhat arbitraary. 

In principle, any model which adequately describes the flow 
field in porous media can be used. This possibility of multiple 
selection, however, , is limited by practical considerations that 
only relatively simple geometrical entities should be used. 

The models employed various investigators can be grouped 
into three categories: capillaric cx)llector:: (Payatakes et.al. ^ 
1974a; Hung and Tien, 1976), spherical collector (Yao et.al. , 
1971? Spielmen and Pitzpatric, 1973? Payatakes et.al. , 1974a, b; 
Rajagopalan and Tien, 1976, 1977), and the constricted tube 
collector (Payatakes £t,^, , 1974c, <3) , 

Using the collector concept, the retaition of particulate 
matter within a filter medium can be considered in terms of 
the deposition of particles from a suspension flowing past 
collectors of specified geometry. The mechanisms of deposition 
are known to be inertial in^action, sedimentation, interception. 
Brownian diffusion, and straining. For liqyid systems, the 
effect of inertial impaction is negligible, and the relative 
importance of other mechanisms depends upon a number of 
variablies, the most significant of which is perhaps particle 
size. Although explicit rate expressions of deposition for 
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some of the individual mec±ianisms are available, in most 
practical situations deposition takes place according to more 
than one mechanism. 

The most versatile method for estimating the rate of 
particle deposition is the so-called trajectory calculation. 

As particles move toward a collector, their trajectories deviate 
from stream lines, and some of them may intersect with the 
collector. If one assumes that deposition occurs once a 
particle makes contact with the collector, the particle depo- 
sition flux can be estimated if particle trajectories are known. 
These trajectories can be determined in turn from the appro- 
priate equations of particle motion, with the knowledge of the 
forces acting on the particle. 

The idea of determining the rate of particle deposition 
from particle trajectories was first advanced 50 years ago in . 
connection with air filtration (Sell, 1931; Albrecht, 1931). 
However, the possibility of extending this concept to liquid 
filtration was recognized only rather recently (O'Melia and 
Stumm, 1967) , 

A detailed discourse on the determination of particle 
trajectories is beyond the scope of the present thesis. To 
initiate trajectory calculation, the following procedures have 
to be followed: 

1, Specification of collector geometry, collector size and 
size distribution* 

2, Specification of flow field around or within the specified 
collector, and 
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3, Specification of forces acting on particles under consi= 
deration, 

Thus/ a variety of oDmbinations of conditions cai be used 
for tile calculation of collection efficiency# a fact which 
accounts for the differences among the various studies on this 
subject carried out during the past decade, A summary of these 
studies is given in Table 2,2, 

The particle trajectories are obtained from the integration 
of the appropriate equations of particle motion# which are 
obtained from the balance of forces and torques acting on the 
particle. The exact fonn of these equations varies with the 
geometry of the collector# but they are first order differential 
equations because of the omission of the inertial force. 
Furthermore# the Brownian diffusion force cannot be included 
owing to its stochastic nature. The forces which have been 
considered are gravity, hydrodynamic drag (with or without wall 
correction)# and different kinds of surface forces. The effect 
of interception is included as a boundary condition to determine 
particle deposition, 

Ihe pertinent equations of particle motion# irrespective 
of their specific forms# contain a large number of dimensionless 
parameters# a listing of whidh is given in Table 2.3. In 
general# these equations do not permit analytical solutions but 
can be solved effectively with a number of algorithms such as 
the fourth-order Runga-Kutta method or the Adams-Moulton method. 
However# for two cases# isolated sphere model ( Rajagopalan and 
Tien# 1977) and the sphere-in- cell model (Rajagopalan and 



Tabiie 2,2. Summary of results on trajectory calculations for deep bed filtration 
(Adopted from Tien and Payatakes, 1979) 
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Table 2,3. Dimensionless groups appearing in trajectory analysis 


Name 

Symbol 

Definition 

Double layer group 

^DL 

^P 

Electrokinetic group 

No. 1 



Electrokinetic group 

No. 2 

^E2 ■ 

2 ’J'Oi ^02'^^'^ 01 ^ ’^02^ 

Gravitational group 


2ap ( Pp - p )g/9ii Vo, 

London group 

^^Lo 

H/9 11 11 a^ V„ 

in 

Relative size group 

Nr 

a /a 

P c 

Retardation group 

^Rtd , 


Brownian diffusion group 
(Peclet n ’nber) 

Npe 



* Here, < is the double=layer reciprocal thickness, i is the 
dielectric constant, and are the surface potentials of 

particle and grain respectively, is the interstitial 
velocity of filter bed, H is the Hamaker's constant, a^ is the 
collector radius, > is the waveloigth of electron oscillation, 
and is the particle diffusion coefficient in liquid. 
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Tien^ 1976) , collection efficiency values obtainec3 from numerical 
calculation have been correlated empirically with relevant dimen- 
sionless parameters^ thus overcx^ming the major objection about 
the use of trajectory calculation for the estimation of collector 
efficiency • These expressions are 


Sphere- in- cell model (Rajagopalan and Tien, 1976) 

’lo = °-’'2 


for < 0,09 


Here, 


A. = 2( 1 - P ) /w 


w = 


( 1 F^) 

2 =. 3P + 3P^ - 2P^ 


Isolated sphere model (Rajagopalan and Tien, 1977) 




( 2,18) 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


( 2 . 22 ) 


The definitions of the dimensionless groups are given in 
Table 2,3, These expressions are based on results obtained 
under the assumption that the surface interactions ( van der 
Wall’s force and double layer force combined) are favourable, 
Ihe collection efficiency obtained from trajectory 
calculation is often referred to as the initial collection 
efficiency or dean collector efficiency. Significant partide 
deposition may diange the collector configuration, alter its 
surface characteristics, and cause the flow field around the 
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collector to be different from that around a clean collector to 
such an extent as to render the collection efficiency calculated 
with the omission of these factors to bo erroneous, Ihe 
principle of trajectory calculation, however, remains valid. 

The problem is that of taking into account all these changes in 
the formulation of the trajectory equation. 

At this stage it is important to note that Eq. (2.18) 
based on trajectory calculations, does not include several 
mechanisms which may be operative ■ in deep bed filtration. 

However, these mechanisms are incorporated in computing the 
transport efficiency of particle collection either empirically 
or based on theoretical calculations (Hall, 1957; Maroudas and 
Esienklam, 1965; Agrawal, 1966; Yao, 1968; Craft, 1969; Cookson, 
1970; Fitzpatrick, 1972; Clint et.al., 1973; Wenk et.al. , 1977), 

2.2,2 The Attachment Step 

The attachment step in filtration is studied using 
the similarity between coagulation and filtration. The resemb- 
lance between coagulation and filtration can be justified from 
the findings of Habibian ( 1971) . Both process can be envisicmed 
as involving a " transport step and an attachment step. 

The overall efficiency of coagulation can be determined 
from the rate at which particle collisions occur multiplied by 
a "collision efficiency factor" which represents that fraction 
of the total number of collisions which produce inter-particle 
attachment. It is plausible that the overall removal efficd^cy 
of a filter may be represented by the rate of transport of 
suspended partides to the filter grains multiplied by a 
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"holding efficiency" which represents the effective portion; 
of contacts, 

The attachment st^ is a colloidal-chemical process# and 
unfortunately the present knowledge of colloid chemistry does 
not provide a quantitative model of this step. Therefore# the 
similarity between the attachment step in filtration and 
coagulation cannot be shown in mathematical terms. The concept# 
however# becomes more evident by considering that in both 
processes the particle should be made "sticky" so that they 
may attach to each other or to the filter grains upon contact. 
This means that the potential energy barrier between the parti- 
cles or between the particles and filter grains should be 
reduced or nullified. 

At. the beginning of a filtration cycle the energy barrier 
between the particles and filter grains can control the attach- 
ment of particles to filter grains. However# as filtration 
proceeds# the suspended particles coat the filter grains and 
further removal occurs primarily the attachment of particles 
in suspension to particles already deposited on filter grains. 
Therefore# shortly after the commencement of a filtration 
cycle# the magnitude of the energy barrier between particles 
in suspension and particles deposited on filter grains controls 
the particle attachment in the filtration process. This same 
potential also governs interpartide attachment in the coagu- 
lation process. The role of attachment step in filtration 
has jbeen c|ualitatively studied (Habibian and O'Melia# 1975) . 
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However^ the quantification of this step remains out of reach 
with the presQit state of knot-Jledge. 

2.2,3 Dynamics of Deexa Bed Filtration Based on 
Theoretical Approach 

The retention of particulate and colloidal matters 
within a filter bed causes changes in the microstructure of the 
filter medium, As a result, the behavior of the flow of fluid 
through the medium and the deposition of particles from suspen-=- 
sions to filter grains can be significantly affected. 

Qualitatively speaking, the major changes resulting 

from particle deposition involves the effective size and 
geometry of filter grains, the surface characteristics of the 
grains, the porosity, and the effective porosity of the bed. 

The last quantity is defined as the fraction of the connected 
void space which is open to suspension flow. This is an impor- 
tant variable in the understanding of the effect of deposition 
on filter performance, since particle deposition may result in 
the blocking of certain flow passages of the bed, A complete 
filtration theory should provide information on these changes 
as functions of operating and system variables. Once these 
changes are known, it would be possible at least in principle, 
to account for these changes in predicting filtration rates 
and pressure drop requironoits. 

The development of modem filtration theories, of course, 
has not reached such a stage as to permit a complete and 
quantitative description of the effect of particle deposition. 
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A number of investigators (Camp, 1964; Mohanka, 1969; Deb, 
1969; Ives, 1969; Sakthivadivel et,^, , 1972) have presented 
relationships between the local pressure gradient of a clogged 
filter bed and the extent of particle deposition. The starting 
point was usually the carman-Kozeny equation. The effect of 
deposition was assumed to increase the grain size and to 
decrease the local porosity on the basis that deposition over 
filter grains is uniform, resulting in the formation of a 
smooth deposit layer outside the grains. In general, their 
results can be expressed as 


^p/^z 


M a. 

= T Cl - p, cj) ^ 

i=l 


( 2.23) 


The number of terms involved M, the coefficients pj^* 5 ,and the 
exponents a^^'s depend upon the particular model used to charac- 
terize the filter bed, similarly, the removal efficiency of 
filters as the filtration progresses is related to the specific 
deposits and the porosity of the deposits. 

The selection and application of these expressions for 
predicting filter performance cannot be made without additional 
evidence oic assuraption on deposit morphology . A complete 
solution to the problem requires a theory which provides detailed 
information about the depositicai process, including the forma- 
tion and growth of particle deposits as a function of time. 

Such a theory is not yet available, although some progress has 
been made along this direction. Payatakes and Tien ( 1964) and 
Payatakes ( 1977’! have presented mathematical models for aerosol 
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particle deposition with, dendrite like pattern and related 
collection and pressure drop as functions of the 

degree of deposition, it is likely that the same approach, can 
be applied to hydrosol systems, Tien 1977) and Wang 

et,_^, ( 1977) argued that d^osition phenomena, in geaieral, can 
be studied in terms of the two properties whici"! are characteri- 
stics of deposited particiesand approaching particles, respec- 
tively, The shadow effect refers to the fact that once a particle 
is captured by a cxjUector, certain parts of the collector 
surface near the deposited particle would no longer be accessible 
to approaching particles C and those areas are called shadow 
areas). The degree of uniformity of deposited matter (or the 
la(dc of it) is a direct consequence of the magnitude of these 
shadow areas. The morphology of particle deposits can be 
determined by the fact that the spatial distribution of particles 
in a dilute suspension follows Poisson’s law and the presence 
of deposited particles allows approaching particles to be 
collected by the collector as well as by the already deposited 
particles,. 

This principle was adopted by Beizaie (1977) in his study 
of deposition on single isolated spherical collectors. His 
results indicate that the deposition process, in general, 
consists of three steps. The first, or the initial stage, may 

f 

be called the dean collector stage dudng which tte presence 
of deposited pajrtides have negligible effect. The second stage, 
or the dendrite growth stage, is characterized by the formation 
and growth of particle dendrites. Daring this stage, an 
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increasing amount of particle deposition is affected fay deposi- 
tion on deposited particles, Tixe third stage of the process is 
termed the open structured solid growth stage. Particle colleo- 
tion is entirely due to deposition on deposited particles, and 
individual particle dendrites become less distinguishable as 
they become merged and interconnected, A direct application of 
Beizaie's results to filter clogging, of course, would be 
inappropriate, since the use of the isolated sphere model 
precludes consideration of the effect duo to the presence of 
neighbouring grains. However, with the use of a more realistic 
porous media model Beizaie* s method on deposition dynamics can 
be applied advantageously to the study of filter clogging. 

With the same concept of doidrite formation, Ali ( 1977) 
proposed a formulation for the change in single collector 
removal efficiency as follows 

d_ 2 

T) = T)a + N ri a (2.24) 

r p p dj. 

Here, is the single collector removal efficiency at any time 
after the start of filtration process, N is the number of retained 
particles acting as collectors, ti is the collision efficiency of 
the clean media grain, a is the ratio of particles adhering to 
particles colloiding with the filter grain, d^^. is the diameter 
of the media grain, and are the corresponding -n, a values 

for particle to particle collision and attachment respectively. 

Further, for estimating the value of N at any time, Ali 
(1977) proposed the following equation: 


8N 

3t 


P” ^ a V I 


(2.25) 
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Here, p”is the fraction of the particles renoved which will act 
as collectors, n is the local cxincentratjx)n of particles in 
suspension, and is the superficial velocity. 

Since, n is a function of toth filter depth arid time, N will 
also be a function of both filter depth and time. The above 
expression fails to consider the variation of N along the filter 
depth. Further, Eq, (2.25) indicates that N would keep on 
acctimulating infinitely and so also the t]^, and hence the overall 
efficiency of the filter bed, which is very unlikely. This 
descripancy is because of neglecting the effect of shadow areas 
in formulating the aforementioned model. 

2,3 Summary 

Recent advances in the study of deep bed filtration have 
shown that inspite of the inherent complexities of the process, 
the problem can be studied fruitfully based on fundamental 
considerations. Resurgent research efforts of -the past two 
decades have already yielded important results of practical 
significance. The review presented, examines both the phenomeno- 
logical and theoretical studies conducted during recent years. 

It reveals that deep bed filtration is a challenging and 
fascinating subject of study and warrants further attention and 
effort frc»i\ scioitists and engineers. 

The main problems which are relevant to a fundamental 
understanding of deep bed filtration are the nature and the 
conditions leading to the retention of particles throughout a 
filter bed, the ciiange of the filter media structure due to 
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deposition, and its effect on filter performance, h reasonably 
complete understanding of the pertinent phenomena is essential 
for the establishment of a comprehmsive deep bed filtration 
theory which can be used as a basis of rational design. 

It appears that the deposition phenomena, in general, can 
be studied in terms of the two properties which are character- 
istic of deposited particles and approaching particles. The 
morphology of particle deposits can be determined by the fact 
that the presence of deposited particles allows approaching 
particles to bo collected by the collector as well as already 
deposited particles. Thus the formulation based upon this 
concept should be able to predict the filter performance in 
accordance with the observation made on deep bed filtration. 

It is important to note that such a formulation must incorporate 
the effect of shadow areas. 



30 


3 . SCCPE OF THE INVESTIGATION 

Recent advances in the study of deep bed filtration have 
shown that in spite of the inherent complexities of the process, 
the problot can be studied fruitfully based on fundamental 
consideration. Many attempts have been made towerds the 
simulation of the dynamic behavior of deep bed filter. The 
major limitation of the existing simulation procedures is the 
requirement of extensive pilot plant studies for the determina- 
tion of the certain coefficients such as filter coefficient, 
specific deposits, porosity of the filter be3, self porosity 
of the deposited material etc. Further, the formulations used 
in such simulation techniques neglect the physical and chemical 
characteristics of the particles forming the deposited material 
in the filter. Though the importance of physical and chemical 
properties of tiie deposited matter was realized long back, 
these could not be incorporated in the formulations because of 
the inherent limitation of estimating such properties. In 
particular, the lack of data on particle size distribution has 
hindered the theoretical studies of filtration process and 
prevented the use of suspended particle size as a basic para- 
meter in design. However, the advent of particle counters and 
their successful use in the control of water treatment process 
is expected to overcome this major difficult^f (Tate and Trussell, 
1976; Wright, 1979). Consequently the present investigation 
is undertaken with the basic motive of developing a conceptual 
formulation for predicting the filter performance. 
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The work is concerned with the development of simulation 
procedure capable of predicting the dynamic behaviour of a 
deep bed filter over the broadest practicable range of the 
pertinent variables involved^ including^ in particajilar/ the 
suspended particle size The method presented here is developed 
by assimilating available quantitative results relating to 
porous media flows* and to filtration* in the cont^t of conce- 
ptual framework. The main hypothesis can be stated as follows; 

"The removal of suspended particles in filtration is a 
two step process. The first step accounts for the transport 
of su pended particles to the filter media grain or already 
deposited particles. The second stop incorporates the 
attachment of such particles to media grain or previously 
deposited particles. This assumes that the deposited parti- 
cles or the retained particles may also act as collectors." 

This picture of the process is supportable on the basis of 

observed filter behaviour and does* moreover* provide the means 

for integrating results into the overall simulation scheme, in 

general the objective of this research is to formulate and 

verify the new conceptual model with aforOTentioned hypothesis. 

Specifically* the present study is directed to the following 

aspects of the overall problem. 

1, To study the role of suspended particles retained in a 
f ihter bed with respect to ronoval efficiency and head 
loss development* and to investigate the influence of 
the size of these particles on filter performance. 

To propose a mathematical model for removal efficiency 
following practice in aerosol and water filtration 
using collector theory. 


2 
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3. To modify the Carman-Kozen y head loss equation for 
clean filters to predict the head loss development as 
the filtration progresses. 

4 . To verify the proposed formulation for removal efficiency 
and head loss developed using existing data and the data 
obtained in the present investigation , and thereby 
evaluating various physical and chemical parameters 
deciding water and wastewater filtration as an unit 
operation. 

5. To simulate the filter performance for various physical 
and chemical parameters based on proposed mathematical 
model . 

6. To propose a mathematical equation for the depth of 
clogging front which would give the depth of floe 
penetration relating to head loss developed as the 
filtration progresses. 

To study the temporal variation of the depth of clogging 
front for various physical and chaaical parameters 
influencing filtrartion process. 


7 . 
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4. MODEL FORMULATION 


Filtration is a complex process due to the diversity of 
the mechanisms involved, A reliable mathematical description of 
the removal of suspended particles within a filter bed is nece- 
ssary to obtain a reasonable understanding of the mechanisms 
involved. The prerequisite to such an approach is the knowledge 
of the mode of deposition of suspended particles within a filter. 

The mode of deposition of suspended particles within a 
filter bed is believed to be in the form of chainlike particle 
dendrites, Tien et,^, ( 1977) have presented a theory for the 
formation and growth of such particle dendrites on a collector 
pltaced in an aerosol or hydrosol system. They attribute such a 
formation to the finiteness of the particle size which creates 
shadowing effect. Due to the formation of a shadow area, part- 
icles find it difficult to deposit in this shadow area and 
hence tend to attach to the previously retained particles 
rattier than the filter grain. Yao et.al . ( 1971) have also 
mentioned the possibility of such particle attachment in water 
filtration. However, the possibility of the existence of long 
particle dendrites in water filtration is a subject of contro- 
versy among the researchers involved in modelling the deep bed 
filtration. Nevertheless, the fact that the filtration effi- 
cienqy changes as the filtration progresses clearly indicates 
that the deposited pairticles play a significant role. 

Most. of the existing mathematical models for predicting 
the dynamics of deep bed filtration attribute changes in filter 



34 


performance to the changes in bed porosity and the self 
porosity of the deposits. These models neglects the physical 
and chemical characteristics of the particles fojcming the 
deposits in the filter. The only term v 7 hich may represent the 
characteristics of the particles is the self porosity of depo- 
sits, For equal amounts of removal, however, it is intuitive 
that self porosity of the deposit is independent of particle 
size. In other words, the head loss resulting from equal removal 
of two suspensions of different particle size will be the same, 
whereas the experimental data of Habibian (1971) show the reverse. 
The filtration of suspension having particle size of 0,1 pm 
resulted in much higher head loss than those having 1 and 7,6 pm 
particle size. Further, if decreased porosity and pore size are 
the main important factors responsible for the changes in 
filtration efficiency, then filtrate equality should continuously 
improve, provided there is no reentrainment of the deposits 
particles. However, many experiments conducted under different 
set of conditions by Habibian (1971) show the deterioration in 
effluent equality right from the beginning of the filtration run 
when the reentrainment of the deposited particles is not 
expected. This discrepancy between experimental observations and 
predicted filter performance using existing models is due to 
the inability of these models to incorporate the physical and 
chemical characteristics of the particles forming the suspension 
into the equations for particle removal and head loss. 

In the present research a new conceptual model has been 
formulated for predicting the filter effluent quality and head 
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loss development as the filtration progresses relating to 
various physical and chemical parameters. The model is formu- 
lated with the concept of deposited particles acting as collec- 
tors taking into consideration the effect of shadow areas. The 
model formulation is described in the subsequent sections. 

4,1 Terms Used 

For the purpose of clarity, some of the terms are described 
prior to the model formulation. 

Single Collector Efficiency: In the filtration literature, 

the term single collector efficiency is used interchangeably to 
represent the single collector removal efficiency without 
distinction. Most of the previous work relating to the use of 
this term deals with the transport aspect of the filtration 
process. Therefore, to be consistent, the term single collector 
efficiency represents the single collector collision efficiency. 

Single Collector Removal Efficiency; The term single collector 
removal efficiency is used to irepresent the actual removal by 
the collector accounting for both the transport and its att^clanenrt 
steps. 

Single Collector: A single collector herein is referred to 

the assembly consisting of a filter grain and that number of 
particles attached to it which act as collectors. At the start 
of the filtration process ( i.e., at time equal tx) zero when no 
pairticles have been attached to the filter grains) the assembly 
consists of only the filter grains themselves. Ihe term single 
collector used herein is, therefore, inconsistent with the use 
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of the term by other investigators who use the single collector 
to represent a filter grain in studying the clean filter 
performance. 

4.2 Assumptions and Their Justification 

To study the role of retained particles in a filter^ a 
mathematical model for the filtration of dilute suspensions 
through a filter bed is developed for a monosized filter media 
and particles / botti spherical in shape. To facilitate mathe- 
matical manipulation/ following assumptions are made. 

1, The flow around the collector remains unaltered as the 
particles are removed. 

2, The diameter of the single collector is equal to the 
diameter of the filter grain and that the diameter of 
the retained particle acting as collector is equal to 
its own size. The assumptions are justified as long 
as the filter media is much greater in diameter than 
the diameter of the suspended particles (Yao, 1968). 

3, The porosity of the filter bed is considered to remain 
fairly constant as the filtration proceeds. This 
assumption is generally made by other investigators in 
water filtration in dealing with the removal efficiency 
of the process. Herzig et.al . ( 1970) / discussing the 
assumption/ have stated that it may result in an error 
of the order of 0.3 percent in dilute suspensions. 
Pavatakes et.al. ( 1974b) /Who used capillaric and 
Brinkman^ s model to calculate the head loss increase 
during the cx 3 urse of filtration/ have reported that the 
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calcxilated head loss values were less by two to three 
orders of magnitude then the observed ones even taking 
into account the decrease in the bed porosity. This 
indicates that in dealing with dilute suspensions, the 
assumption of constant porosity is not critical, at 
least for some time after the start of filtration. 

4. The flow through the filter is laminar. Justifying 

this assumption in water filtration, Yao ( 1968) states 
that for filtration velocities from 0,136 to 0.272 cm/sec 
and media sizes upbo 1.0 mm, the Reynolds number based 
on media size varies approximately from 0.5 to 5, 
whereas the transition from laminar to turbulent may 
begin at Reynolds number from 10 to 113. 

4,3 Model Representation of Filter Bed 

To provide a physically realistic model as the basis for 
quantitative descriptions of the various processes of interest 
which occur in the porous medium, the concept of Unit Bed 
Element (UBE) and the unit collector proposed by Payatakes et . 
al. ( 1973) is used. Thus, the filter bed is postulated to 
consist of a number of unit bed elements, with a physical height 
of AL, connected in series. Each UBE in turn is assumed to 
consist of a number of particle collectors . The capacity of the 
unit bed to remove particulate matter can thus be defined as the 
fraction of suspended particles entering the UBE removed. The 
general idea of depicting a filter bed as a series of ube's is 
illustrated in Pig. 4.I. 
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Fig. 4.1. Representation of filter bed as an assembly 
of collectors. 
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Specification of geometrical configurations for the 
collectors present in a UBE is admittedly somewhat arbitrary. 

In principle, any model which adequately describes the flow 
field in the porous media can be used. This possibility of 
multiple selection, however, is limited by practical consider- 
ations that only relatively sinple geometrical entities should 
be used. The models employed by various investigators are 
reported in section 2.2.1. Because use of the spherical coll- 
ector leads to significant simplifications, it ramains as an 
attractive collector configuration, affoarding considerable 
computational convenience as well as leading to physically 
meaningful results when used in appropriate context. Further 
Rajagopalan and Tien ( 1976) reported that the calculated 
collection efficienqy values based on different collector confi 
gurations turn out to be comparable in magnitude. As a conse- 
quence of the foregoing considerations, it is proposed to use 
the Happel ‘ s porous media model to compute the collision 
efficiency in this study. 

4.4 Hypothesis of Deoosition Process 

The importance of the knowledge of deposit morphology in 
understanding the effect of particle deposition is rather 
obvious. Therefore any theory presented for the performance of 
the process must account for the manner in vdiich the particles 
are deposited within a filter. 

In the present work, the deposition phenomenon is studied 
in terms of the two properties which are characteristics of 



deposited peirticles and appiroachinQ panticles respectively. 

It is postulated that scxne of the retained particles within a 
filter can act as collectors along with the filter grains. 
Further / once a particle is captured by a collector/ certain 
p3tts of the collector surface would no longer be accessible to 
approaching particles. In other words the deposition process 
may be viewed as follows . 

"Initially when the filter bed is clean/ the deposition 
process is controlled by the interaction between the suspension 
particles and the media grain. As the filtration progresses, 
deposited particles cover certain parts of the media grain and 
the deposition is controlled by the interaction between media 
grain and suspension particles, and deposited particles and 
suspension particles for the uncovered and covered portion of 
the media grain respectively. For the particles which are in 
direct contact with the media grain, the transport step would 
be controlled by the size of the grain and suspension particle 
while the attachment step would be controlled by the character- 
istics of deposited particle and approaching particle for the 
interaction betweai deposited particles and suspension parti- 
cles, However, the portion of the retained particles which 
are acting as collectors and are deposited on the retained 
particles, both transport and attachment would be controlled 
by the properties of deposited particles and suspension 
particles." 

4,5 Development of the Model 

4,5,1 Prediction of Rai:\oval Efficiency 

Following the practice in aerosol and water filtra- 
tion ( Friedlander, 1958; Yao et.al. , 1971), the single collector 
efficiency t)^ can be obtained as follows 

_ Rs-te at which pa articles _strike_ i^e oDl lector 
^c Rate'”a't whidi particles move toward the collector 


Considering a spherical collector of size *d^', the projected 
area noamaal to the flow is 


2 
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The rate at which particles move toward the collector is, 
therefore 


= V n 
o 


% ,2 

T Q 
4 c 


Since the collector consists of a filter grain and that number 
of particles attached to it whicih act as collectors (N), the 
single cxjllector efficiency can be rewritten as follows. 


Rate at which 
particles strike + 
the filter grain 


rate at which particles 
N X strike a retained particle 
^acting as collector 


V n -7 d 
o 4 


2 

c 


Rate at which 
particles strike 
the filter grain 


V n 
o 


TC 

4 



Rate at which particles 
strike a retained part- 


+ N X (^) X collector 


„ 7t ,2 
V 4 


or 


r\ 

c 


ri + N T)^ 



(4.1) 


Accounting for the attachment step, the single collector removal 
efficiency ‘T)^' can be expressed as 


d 2 

\ = "na N dh/d^ + Nti (^) 

a 2 

+ ( 1 - P ) N T)„ a_ ( (4.2 
•r V * Kc 'p p d^ 

Here, p is the fraction of retained particles acting as coll- 
ectors which would cc«tribute for the surface coverage of media 
grains, and a and a are defined following the practice in 

coagulation as follows. 
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Number of collisions between filter grain and 
^ _ ®-®P®^^^i?.-P5?^tiGles whi^ attaclunent_ 

Number of collisions vfiaicii occur between filter 
grain and particle 

and 



Num^r of collisions between retained particles 
acting as collectors and suspended particles 
.ip a_ttactoent 

Number of collisions whicii occur between retained 
particles acting as collectors and suspended 
particles 


In other words / a and a can be called particle to filter 

Jcr 

grain and particle to particle attachment coefficients respec- 
tively . 

In Eq. (4.2) the factor [( d^ - B N d^)/d^l is introduced 
for the reduction in surface coverage of the media grain which 
would not be available for contact of suspended particles. 

The second term in Eq. (4,2) represents the removal efficiency 
because of the retained particles acting as collectors vhich 
contribute for the surface coverage of the media grain. For 
such retained particles the collision efficiency is t] and the 
contact efficiency is a^. Ihe third term represents the ronoval 
efficiency of the retained particles acting as collectors 
which would not contribate for the surface coverage of the 
media grain or which are not in direct contact with the media 
grain* Ihe collision and contact efficiencies for these 
particles will be and respectively. 

The single collector efficiency of a filter grain is 


expressed as 
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T) = p^^ticles^ strike the filter gr p in 


v n |d2 
O 4 G 


or rate at which the particle strike the filter grain is 


n V n f 
' o 4 c 


If a is the particle to filter grain attachment coefficient, 
the rate of renaoval of particles by a single filter grain is 


r)a 


-*-i- >-t ^ p ^ 

V n — d 

o 4 c 


Similarly, the rate of removal of particles by retained part- 
icles acting as collector is 


r| a V n d^ 
P p o 4 p 


The change ‘in number of collector particles on a single filter 
grain can be written as follows. 


Change in number of collector particles on a single filter grain 
= Rate of removal of particles by single filter grain 

+ N X rate of removal of particles by retained particles 
on filter grain acting as collectors 

-NX rate of parti ciLes removed which wmld contritute 
for the reduction in effective area of the 
filter grain 

2 

- N X rate of particles removed which wculd contribute 
for the reduction in effective area of the 
retained particles acting as collectors 
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n ^2 


r,a v^n 5 N t, v^n f + ( l.p^)N v^n f d^ 


« 'Pc ^ V f ^p + (1-p^) np«p v^n I d^} 


- Ti a V n I d^ P 

P P O 4 p ^P 


11 ^2 
P 


[ria u ^ N tiap I + (i.p^)N f d, 

^ “p V, f + ( 1-p^) ^papV^|d2}p^ 

.2 .. % 7 


N 


TL a„ V d 6 In 
, p p o 4 p f^p -I 


[Cj^ + N c^ + (i - P^lCj )- N Cj + (l-p^)C2}Po 


~ N C2 p ] n 


= [C^ + N(C3 - C3) - N^ C4] n 

Here, p is t±ie fraction of removed particles which would 
P 

contribute for the reduction in effective area of the retained 

particles acting as collectors. 

^2 

Cl = n d v^,. 5S , Cj = rip dp Vp f d" ; 

=3 = ^ °2 ^ '■ '^4 = '^2 Pp • 

p 

Considering an extronely tbin UBE of thickness ‘AL' of 
the filter bed, and making the mass balance for N, the following 
expression can be c^tained» 



45 


aN Sn r , o 

at aiT “ + n(C3 » ^^03) - n^c^] n (4.3) 

The removal efficiency of the filter bed can be determined by 
considering an UBE of thickness 'aL' of the filter bed. The 
number of single collectors in this layer 

A( 1 - F) A L 
c 

^ ~ 

Here^ A is the area of the filter bed. 

The removal by a single collector of removal efficiency t) 
from definition 


Therefore the ran oval by the UBE of thidcness 'AL‘ 


3 

2 




A V n 
o 



Making the mass balance for the suspaision as it flows 
through the UBE under consideration, a mass balance equation 
can be written as 


Rate of accumulation = Rate of input - Rate of output - Removal 

If 'n‘ is the number of particles in suspension entering 
the UBE, the changes in number as they flow thrcxigh the UBE is 

3 n . _ 

" + sT 


and 
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Rate of input = A v n 

o 

Rate of output = + 

o 9 L ^ 

Rate of accumulation = a al — 


In matiiema tical feinas the mass balance ec3uation becomes 


A AL 


= A V^n - A V^(n + — AL) - I a V n i 


( 1-F) 


3n 

at 

or 



j. V in ^ ^ 1-”^) 
3t^ oaL^2 d 


V n 
o 


o 


(4.4) 


Eqs, (4.2), (4.3) and (4.4) comprise a mathematical model 
to describe the removal efficiency of a filter bed in time and 
space. Their solution results in the removal efficiency of a 
filter bed. The exact solution of the above equations is not 
possible. Simplifying techniques are used to obtain a reason- 
ably accurate solution. In this r| and n are considered as 

r 

step functions rather than continuous functions of time and 
depth. A new algorithm has been developed to solve these 
equations employing numerical methods. The details of the 
solution technique, are presented in Appendix A. 
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4,5,2 Prediction of Head Loss Development 

The development of head loss (h^) is based on 
the well established Carman-Kozeny equation for head loss in a 
clsan filter ( and can be written as 



„V(1=F)^ A 2 

K' E ^ „ p ,_cv 

P g „3 ''V ^ 

Fq C 


(4.5) 


Here^ and are the surface area and volume of filter 

grains, and K* is an empirical constant. 

As the filtration proceeds, the suspension particles are 

deposited in the filter bed and thereby increases the surface 

area and volume within the filter. If the retained particles 

in time t contribute additional surface area A and volume V 

P P 

to the filter, the Carman-Kozeny equation becomes 



H _jS 1 •” F) 


A + A 2 

(-C .E) 

V + V 
C P 


,, V .2 Tl d^ N + P' 7C d^ N 2 

» M _£ 1 ", r_ __ P , c ,, -P. P-] 

P 9 ^ f d^ N + I d^ N 

6 c c 6 P P 


(4.6) 


Here, N and N are the number of filter grains and retained 
c p 

particles in the unit volume of the bed respectively, h^ is 
the head loss after time t, and P' is that portion of the 
particles which contritute to the additional surface area. 
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The Eq, (4,6) can he rewritten as 


36 
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£ 

g 


M d 2 

0 1 + 6 * 

[- =-. 5 - ^3 ] 

c c 


(4.7) 


Due to non—uniform deposition of particles in different layers 
^f 

of the bed/ — in each layer will be different. The filter 
depth may be divided into a reasonable number of layers and 


h, 


h. 


can be calculated in each layer. The sum of r-^ for all 
Li L 

layers will result in the total head loss across the filter 


bed. N in each layer can be calculated from the renoval 
P 

efficiency . 


4.6 Modification of the Model to Incorporate Varying 
Porosity 

The model developed in the previous section, is based on 
the assumption that the porosity remains constant throughout 
the filtration process. Though this assumption is justified 
in case of dilute suspensions, to be more realistic the vari- 
ation in porosity should be considered in the model formulation. 
The variation in porosity is incorporated by assuming that 
Eqs. (4.2), (4.3), (4.4), and (4.7) are valid only for a very 
small time interval At in which the porosity remains cx>nstant. 
After every time interval the porosity value is modified by 

AF . as follows: 

J 

p % =i F — AF . ( 4 ,8) 

j+1 0 3 
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Here, ±S the porosity in the ( j+l)** tin« internal ana 

JFj Is the change in porosity in the (j+1)^ time interval 
compared to initial porosity. 

The change in porosity is computed by the volume of the 
deposited particles per unit volume of the filter bed as 
expressed by the following equation 


AF. 


\ (n. - n ) f d^ 

1 1 j e^ 6 p 


(4.9) 


Here, n. and n are the influent and effluent particle cx>ncen 

j i 

trations respectively . 

•(- 4 ^ 

The porosity for the (j+1) time interval can be written 


as 


F = F 

j+1 o 


-I (n. - n )| 

i ^ ^ ^ 


(4.10) 


Ihus taking the appropriate value of porosity, the assumption 
of constant porosity can be eliminated in computing removal 
efficiency and head loss d velopmoit by Eqs, (4.2) to (4*4 ), 
and Eq. (4,7) respectively. 


4.7 Modification of the Model for Varying Particle Sizes 
in the Suspension 

Ihe model developed in Section 4,5 is valid only for 
monosized particles in the suspension. However, the suspended 
pairticles presesnt in either water or wastewater are of varying 
sizes. In orxjer to ajply the mcdel develop^ to the realistic 
situations certain itodificatioi^are essential, Ihe model for 
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monosized particles in suspension may be extended to varying 
particle sizes as follows. 

Consider that the influent suspension contains m particle 
sizes having diameters d , d , in concentrations 

^ 1 X ^ 


/ . • • / n . Obviously 
1 2 m 


'n 


n + n +.»».+n = n 

'"l ‘='2 Om ° 


or 
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In = n 
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For the j particle size Eq. (4.2) may be written as 
> ■ o m 


' = Vj ‘‘’c - Pc ''j Sp 


^ )/de 

J 


m 


m 


+ [ I N. I T} 


oc^ 


u - 'n L f 4 

i=i j=i J Pi,j 


d 0 

P. 2 


m m 

+ V . “P. 




p^ 2 J-m 


i=l j=l ^i,.j "'i^i c U=1 


(4.11) 


(4.12) 


til* til 

Here^ subscripts i, j refer to the i and j particle sizes 
respectively . 

Ihe equation for computing the number of retained particles 
acting as collectors of the j pap-ticle size may be obtained 
from Eq, (4,3) as follows. 


3N. 3N. 

I gt o 3 L 


C, n. + 

ij i 


m 

x=l 




m 


I C 

i=l 


n 


t m 




(4.13) 
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Here , 


^ j ^ 4^ 




' ^2, . ^c “ ^c^2, . ' ^4. .“ ^2. P 


i/J 'Pi . x,J 


c i. . 

3-/3 


4. . “ ^2. . ^'p 
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Similarly Eq, (4,4) for the J particle size may be written as 


■ '. ""i + V ? + 3 L1._- .F) ^ ^ ^ 

j 9t ^ 0 9L *2 d ^o ^r = ° , 

c -» j h=i 


( 4 . 14 ) 


Eqs. (4,12), (4,13), and (4,14) comprise a mathematical 

model to describe the removal efficiency of a filter bed in 

tin 

time and space for the j particle. The overall effluent 
quality may be expressed as 


n = n. + n^ + , . . + n„ 

X X rn 


n = in. 


(4.15) 


The head loss development may be predicted the 
following equation which can be obtained from Eq, (4,7), 

m l^p. ^p- 2 

h , .2 ^ + P' i N‘ '■ ^ d 2 

^ = 36 B , V (4.15) 

L 0 Cf o ^ ^ ct Q 

^ d m p/ V* ^ 

c ... y _ 3, (..1.) 

j=l c c 

The verification of the proposed model for removal 
efficiency and head loss development using the data available in 


the literature and the experimental investigatior> ^^^uct^ 

this research is descrihed in Chapter 6, ^ ^ 


this research is descrihed in Chapter 6, 
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5. EXPERIMENTAL MEIHODOLOGY 


5.1 S cope 

The validity of the model and the proposal based upon 
it were evaluated by conducting appropriate experiments to 
the data available in literature. These experi- 
ments are discussed in detail in Chapter Six of this thesis. 
The effect of different particle sizes in the suspension 
on filter performance was investigated by conducting labo- 
ratory scale filtration experiments in both regions of 
filtration and in the critical suspended particle size range. 
Alum was used for precoating the media and destabilizing the 
suspended particles. Coagulation experiments were conducted 
to estimate the appropriate dose of alum for destabilizing 
the suspended particles in filtration experiments. 

To facilitate the discussion of the apparatus# 
materials# and procedures# these experiments are divided 
into two groups: (1) the filtraticai experiments which were 
designed to investigate the influence of suspendel particle 
sizes and their distribution on removal efficiency and head 
loss development, and ( 2) the destabilization experiments# 
which were conducted to estimate the appropriate alum dose 
to be used in filtration experiments. In subsequent sections 
of this chapter the apparatus# materials# arsi procsiAires used 
in these investigations are described. 
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5.2 Filtration Experiments 

Many factors affect water and wastewater filtration and 
require control in filtration research, Tliese factors can be 
categorized as follows: 

A. Characteristics of the filter bed 

(a) size and shape of the media grains 

( b) bed depth 

{ c) bed porosity and pore structure 

(d) surface characteristics of the filter media 

B. Characteristics of the suspension to be treated 

(a) size and shape of the suspended particles 

( b) density of the suspended particles 

( c) surface characteristics of the suspended particles 
(d) concentration of the suspended particles 
( e) chemical characteristics of the water or wastewater 
(e,g,, pH, ionic strength etc.) 

( f ) chemical characteristics of the coagulant or 

"filter-aid" (e.g., charge, molecular weight etc.) 

C. Characteristics of the operation of the filter 
( a) filtration rate 

( b) direction of flow 
( c) badcwashing procedures. 

In this research^ characteristics of filter bed and its 
operation were kept constant in most of the experiments. 
Emphasis was placed on the characteristics of the particles 
to be filtered. More specifically, the filtration of suspended 
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particles in both regions of filtration as well as in the 
critical range was studied when present separately and in 
different conbinations . 

5.2.1 Characteristics of the Filter Beds and Their 

Operation: 

The filter apparatus (Fig. 5.1) employed in this 
investigation was principally designed to provide for testing 
the proposed model. In few experiments it was decided to 
employ two filter beds having different media to operate in 
parallel. Both the filters received same influent. 

Laboratory filters were prepared from plexiglass tubes 
having an internal diameter of 2.6 on. The use of such a 
small filter was necessary in view of the high cost and 
limited supply of latex particles \^ich were used to prepare 
influent suspensions. It was intended to compare the results 
of the filtration experiments of this research with those of 
Habibian (1971). Therefore, it was decided to use more or 
less the same filter depth and filtration rate that Habibian 
C 1971) employed in his experiments. The filter depth in most 
of the experiments was 15.0 cm. The use of such a shallow 
filter in comparison with the 75 cjti of bed employed commonly 
in a conventional rapid sand filter was necessary to get the 
detectable turbidity in effluent. Moreover, the objective 
of the experimental investigation was to verify the proposed 
mathematical model rather than attributing more significance 
to the extent of removal. 
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> bed porosity of 0.429 and a filtration rate of 0.136 
cju/sec were employed. 'Ibese values are more or less similar 
to those used in conventional rapid sand filters. 

Glass beads (supplied by Li^it and Co,# Azamgarh, India) 
were used as filter media. These beads had mean diameters of 
0.3 mm and 0,5 mm (data supplied by the manufacturers). 

These beads were quite uniform and spherical in shape to the 
extent which can be seen by the microscope ( lox). However, 
in few experiments using clay suspension, sand with Geometric Mean 
(GM) diameter of 0,297 mm was employed. 

The ratio of the diameter of the filter column to the 
grain diameter was greater than 50 and 85 for 0.5 and 0.3 mm 
glass beads respectively. It has been reported that wall 
effects are negligible for values of this ratio greater than 
50 ( Rose, 1951) . 

The filters were mounted on wooden frame. A fine brass 
wire mesh was inserted in the filter bottom to prevent the 
escape of filter grains. Provision for measuring the head 
loss across each filter at various depths was made. Polyethy- 
lene tube manometers, 0.6 cm in diameter, were used to measure 
the head loss. Connections were made with rubber tubing. 

Two manometer lines were ccainected to dimatrically opposite 
taps at eadi level in the filter column opposite taps. This 
avoided erroneous readings of pressure heads due to air 
blocks clogging etc. Normally the two mancOTeters gave equal 
or close readings and in the latter instance a mean value was 
used. A sampling point at the effluent side of each filter 
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was provided. A flow meter (Maximum Flow Rate: 84.6 cc/min; 
Model: 5 GC; supplied by iepl, Hyderabad, India) was installed 
on each sampling tube with an accon^anying screw clamp to 
permit the adjustment of the flow. Test tubes, about 20 cc 
in volume, were used for sample collection. 

5.2.2 Preparation of the Filter Bed 

Elaborate preparation of the filter beds was 
required , Every time the glass beads were first soaked in 
0,1 N NaOH followed by several washings with distilled water. 

To alter the zeta potential, filter media was coated with 
alum as per the procedure described by Sriramulu ( 1975). 

About 275 gm of filter media was taken in a one litre plastic 
bottle containing 550 cc of 50 mg/lit alum solution and kept 
on rotating shaker (supplied by Gansons Pvt, Ltd,, Bombay, 
India) at 40 RPM for 1 hr. At this high dose, alum formed 
micro floes and coated the media. The media was then care- 
fully transferred into the filter column without agitation. 

This procedure ensured a uniform coating of filter media. 

Mater from the clean water overhead tank was then 
passed through the coated filter beds to permit establishment 
of the desired rate of flow through the system. This tank 
was made of plexiglass and equipped with an ovejrflow and 
two outlets, each connected to one of the filters. 

-3 

The distilled water containing sodium bicarbonate ( 10 

-3 

M) and sodium chloride ( 10 M) was employed to establish the 
flow. Thns this water resembled the actual influent to the 
filter with respect to solids concentration, pH and alkalinity. 
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A small centrifugal pump was used to lift this water to the 
overhead tank. . 

After setting the flow, the filter beds were tapped 

slightly to produce the required filter depth and thus the 

required bed porosity. Before introducing the suspension into 

the filters, the "equality" of filter beds so prepared was 

Ah 

checked. This was done by the comparison of the ratio * 

Here, ^h^ is the initial head loss across the filter (cm) and 

Q is the rate of flow thrcugh the filter ( cc/min) . The ratio 
A h 

“ “q^ ^ had values of 0,142 and 0,188 for 0,5 mm glass beads at 

15.0 and 20,0 cn\ respectively. The corresponding values for 

0,3 mm glass bsiads were 0,393 and 0.524, If the measured value 
Ah 

of for any filter deviated more than +0,005 from the 

above values, the preparation of that bed was repeated until 
this criterion was satisfied, in all experiments, the head 
over the media surface vjas more or less kept ccaistant at 

1,73 + 0,01 m during the entire length of the filter run. 

\ 

5.2,3 Characteristics and Preparation of the 
Suspension 

After completing the preparation of the filter 
bed and the setting of the flow rate the suspension to be 
treated could be introduced to the filters. Latex particles 
(Dow Chemical Co.) with mean diameters of 0,109, 1,1 and 7.6 
micron were utilized to make different suspensions. Like 
most of the colloidal particles fcund in natural water and 
wastewater, these latex particles are negatively charged and 
thus provide a useful model cx)lloidal suspension. Furthermore, 



S9 


they are very uniform in size (Table 5.1). This uniformity 
is quite useful in filtration research, because the efficiency 
of packed bed is quite sensitive to the size of the suspended 
particles. 


"^^.ble 5.1, Some properties of latex Particles 


^ rv ^ ^ A. « 

f 

Mean Size, ] 
micron 

1 

Standard Deviation, 
micron 

I No . of 

1 Measurements 

1 

f 

{ Density, 

I gm/cc 

1 

0.109 

0.0027 

318 

1.05 

1.1 

0.0059 ■ 

106 

1.05 

7.6 

2.3 


1.05 


* Data supplied by the Dow Chemical Co. 

The influent suspension was prepared by proper dilution 
of concentrated latex suspension. Dilution was accomplished 
with a distilled water-sodium bicarbonate-sodium chloride 
mixture identical to that employed for setting the flow. 

After equilibration with the laboratory atmosphere the suspen- 

-3 

sions had a concentration of 10 M of sodium chloride and 
sodium bicarbonate and maintained the pH at about 8,5 with an 
ionic strength similar to that of a typical natural water. 

The suspensions of these particles were stored in a 
70 litre polyethylene tank prior to use. The contents of 
this tank were mixed a heavy duty stirrer ( Model Q— 1-125 , 
supplied by REMI, Bombay^ India) to assure uniform particle 
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concentration within the entire hulk of the suspension. 

in few experiments, kaolinite suspension prepared in 
the laboratory as per the procedure described by Salome (1980) 
was used. The ionic strength and pH of this suspension was 
same as that of latex suspension, 

5.2.4 Addition of Destabilizing Agent to the Filter 

Influents 

Because of the nonavailability of Indian made 
cationic polyelectrolytes and the problems involved in importing 
such polyelectroly tes , it was decided to use alum as destabi- 
lizing agent, GR grade alum (supplied by S.M., chemicals, 

Bairoda, India) with molecular weight 630,39 and chemical 
formula Al 2 ( SO^) ^.IBH^O was used. The desired dose of alum 
was added from the stock alum solution (0,5 percent) prior 
to pumping the influent suspension, 

5.2.5 Measuring the Filter Efficiency 

The filter efficiency can be determined by 
measuring the concentration of suspended particles in the 
filter influent and effluent, often the nature and extent of 
experimentation dictate the method which should be employed 
for measuring the process variables. 

A light- scattering method seemed appropriate. For 
convenience, optical density was used as an indicator of the 
amount of light-scattered; a Backman Model IU2 ultravoilet 
spectrophotometer was employed. For kaolinite suspensions, 
optical density was measured ty spectrophotometer ( supplied 
by Systrcaiics, Ahemdabad, India) at wavelength of 440 mjA. 
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The optical density of latex suspensions with particle sizes 
0,109, 1.1 and 7,6 micron were measured at wavelengths of 
270, 490 and 790 mjj,. Calibration curves were prepared by 
measuring the optical density of the known suspensions at 
these wavelengths. The curves proved to be linear for the 
range of concentrations used in this investigation. Accordingly, 
the ratio of filter effluent to influent optical densities 
was used in evaluating the filter efficiency. The optical 
density changes with particle size, thus the existence of 
particles in the filter effluent which might have grown in 
size during the course of filtration can reduce the accuracy 
of the evaluation. However, none of the many microscopic 
observations made by Habibian ( 1971) on filter effluents for 
1,1 and 7.6 micron latex particles revealed the existence of 
any particle aggregates. Thus the reported ratio can be 
considered as the real ratio of particle concentration in 
filter effluent and influent, 

i#ien using 0,109 micron latex, the particles could not 
be observed using an ordinary optical microscope either in 
the filter influent or effluent. In view of the very rapid 
transport of these particles by diffusion, the occurrence of 
particles aggregate in the filter effluent is theoretically 
possible. Since such occurrence could increase the optical 
density (Habibian, 1971), the reported efficiency for 0,109 
micron particles can be slightly underestimated . 
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5,3 Destabilization Experiments 

Destabilization experiments were conducted to determine 
the alum dose to be used in filtration experiments. The idea 
of using alum was to reduce the surface charge of the particles 
so as to enhance the particle=to-particle and particle to 
filter media grain attachment. The procedure used for deter- 
mining the desired alum dose is described in the following 
section . 

5.3,1 Jar Test 

Jar tests were conducted to determine the 

appropriate alum dose necessary for destabilization of latex 

and clay particle suspensions. These suspensions were 

prepared by proper dilution of the concentrated stock 

suspensions with distilled water. Sufficient sodium chloride 

and sodium bicarbonate were added to the distilled water to 

-3 

provide a concentration of 10 M of each of these substances. 
The jar test apparatus (supplied by PHIPPS & BIRD Inc., USA) 
was used for the destabilization experiments. Equal volumes 
of the suspension to be destabilized were added to each of 
the six jars used. Appropriate amounts of the alum was then 
introduced into each jar, while the contents of the jars were 
rapidly agitated. Ihe rapid stirring at lOO RPM continued 
for 2 minutes followed by a period of slow stirring at 30 fpM 
for 20 minutes. The residual light absorption by the 
suspension after 30 minutes of settling was used as the 
index of the extent of particle destabilizaticai. ihe maxxmum 
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possible aluiTi c3osg which did not show any significant 
difference between initial (prior to addition of alum) and 
final (after destabilization) light absorption was chosen 
for destabilization. For each particle size an appropriate 
wavelength was selected (Section 5,2.5). 
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6. MOCEL VERIFICATION 

The success of any mathematical formulation for 
predicting the dynamic behaviour of deep bed filters depends 
upon the extent of quantitative agreement between model 
Predictions and the observations made in laboratory and field. 
Therefore it is essential to verify the proposed model using 
laboratory and field data on deep bed filtration. 

Extensive data on filter performance are available with - 
different sets of operating and system variables. However, 
the lack of information on particle size distributicn stands 
as a major obstacle in utilizing such data for verifying the 
conceptually developed formulations for predicting the filter 
performance, very few studies are conducted in filtration 
research controlling the suspoided particle size, Habibian 
(1971) has conducted extensive studies on filtration using 
different particle sizes in the suspaision. Therefore the 
data presented therein have been exhaustively used to verify 
the proposed model. Besides some more experiments are cond- 
ucted in the present research under differ^t set of conditions 
with special reference to varying suspended particle sizes to 
check the validity of the model in more realistic situations. 

The use of proposed model [Eqs, (4,2) to (4,4), and (4.7)] 
to describe filter performance requires evaluation of the 
terms t). Tip, a, ^ Pp » single collector 

removal efficiency (t) a) is calculated from the exqperimental 
data for a 'clean' filter, using Eq. ( 4.4) at t = O as follows . 
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- 3 ( >P) ^ 

n_ _ " 2 'd7‘ • ^ • na 

® ^ ( 6 . 1 ) 

or 

tia = - in (2-)/[| . . L] (6.2) 

O '-‘g 

Here, L is the depth of the filter bed. 

With this initial oondition established, the constants 11 ^ a^/ 
and Pp are computed using non-linear curve fitting tech- 
niques. Morquardt ( BSOLVE REGRESSION ALGORITHM) is used to 
solve for the coefficients in nail ti variable non-linear 
regression equations (Kuester and Mize, 1973). The details 
of this algorithm are presented in Appendix B. The values of 
T) and r)p are calculated using Happel's flow model- { Happel , 1958). 
The pertinent equations are described as follows. 


T1,I 

^G 


4,04 A 


1/3 


2 s 


( Pg - p) 



( 6.3) 

(6.4) 

(6.5) 


and 

n = \ + nj. -t ^ 

Here, rij, and Vg represent the dimensionless single coll- 

ector collision efficiencies when the re^ective transport 
mechanisms are diffusion, interception aid sedimentation, and 
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is the Peel et number given by the following expression. 


V d 
-.9 .... C 

^BM 


< 6 . 7 ) 


and 


^ p 

The cx)niputed values of the T)^ and t) for the dlffexanfc 

set of conditions employed in the verification of the laodel 

are presented in Table 6.1. The computed value of h using 

IP 

Eqs. (6,3) and (6,4) comes cut to be more than 1,0, Since the 
collision efficiency cannot be more than 1.0, t)^ Is taken equal 
to 1,0, The values of a and a are canputed as follows . 

XT 

a = (T)a)/ifi (6.9) 

The value of p* is computed using Eq. (4.7) in a similar way 
to the computation of t) a , p , and p , The model verifi- 

P p o p 

cation is presented in two different sections using 
1. The data available in literature, 
and 2. The data obtained in the present investigation. 

The validity of the model was checked by the two statis- 
tical parameters, the Coefficient of Correlation (R), and the 
Standard Error of Estimate (S). The expressions for R and S 
are given in Appendix C, The low value of S together with the 
value of R approaching the unity is an indication of the gcxDd 
agreement between model predicticxis and experimental ohser— 


vations 
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Table 6,1, Computed cxdlision efficiencies for clean collector 
using Happel's flow model (Happel, 1958) 

Temperature = 300°K; g = 931,0 cm/sec^; p = 1,00 gm/cc 



1. d = 0.038 cm; F = 0.36; V = 0.136 cm/sec; P = 1.05 gm/cc 

^ 00 p 


0.109 

0.2357x10"*^ 

0.1625x10“'^ 

0.7980x10"^ 

0.7997x10'°^ 

1.100 

0.2400x10"^ 

0. 1655x10*"^ 

0.1709x10“^ 

0,3388x10"^ 

7.600 

0.1146x10”^ 

0.7899x10"^ 

0.4711x10"^ 

0.8061x10”^ 

2. d 

c 

= 0.03 cm; F 

0 

= 0.429; V = 0. 

0 

136 cm/sec; p = 1,05 gm/cc 

0.109 

0,2357x10“^ 

0.1676x10”^ 

0,8062x10“^ 

0.8079x10’^ 

1.100 

0.2400x10'^ 

0.1706x10"^ 

0.17 26x10“^ 

0.3457x10”^ 

7.600 

0,1146x10"^ 

0.8146x10“^ 

0.4759x10"^ 

0.8308x10"^ 

3. 

= 0,05 cm; F 

0 

= 0.429; = 0, 

136 cm/sec; p = 1,05 gm/cc 

Jtr 

0.109 

0.2357x10”^ 

0.6032x10"'^ 

0,5735x10"^ 

0.5742x10“^ 

1.100 

0.2400x10“"^ 

0.6143x10“^ 

0.1228x10"^ 

0.1866x10”^ 

7.600 

0.1146x10"^ 

0,2932x10“^ 

0.3386x10”^ 

0.3081x10-^ 
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6.1 


Modal Verification 
Literature 


Using the Data Avail ahle 


in 


4i.s Stated earlier, cnly Habibian (1971) has conducted 
filtration experiments controlling physical and chemical para- 
meters influeicing the filtration process with particular 
reference to the size of the suspended particles, and hence his 
data is extensively used for the verification of the proposed 
model , 


In one experiment, Habibian operated four identical 
filters C d^ = 0,038 cm, f = 0,36, L = 14,0 cm) at a flow rate 
of 0,136 cm/scc (4.8 m /in /hr), filtering a suspension of latex 
spherical particles (n^ = 11,0 mg/lit, d^ = 0,109 |im, = 

1,05 gm/cc). Filter media were coated with cationic polymer 
prior to use, and polymer was added continuously throu^out the 
runs at a dosage selected to maximize particle- to-particlo 
attachment. Coagulation within the pores of the filter bed was 
reported to be insignificant. Experimental results of the four 
runs wore similar; these have been averaged and plotted in 
Fig, 6,1, 

Ihe sin^e collector removal efficiency ( u a ) is calculated 
from the experimental data for a 'clean' filter. For example, 
n/n^ « 0.6 for the filters in Fig. 6.1. Hence, using Eq. (4.4) 
at t = 0, T)a = 1.44298 x 10"^. With this initial condition 
established, the constants 11 ^ Pp/ and p' are determined 

using non-linear curve fitting technic^es. ihe values of 
1) (T/ p / p , and p' are 0,041, 0,89 , 0,13 x 10 and 0,22 

respectively, Ihus tt and ffp can be ccmputed as 0,171 and 0,041 


respectively , 



HEAD LOSS, cm n/n#, percent 


100 
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80 ^ 



dc = 0.038 cm 
Vo 5 0.136 cm /sec 1 

L =14.0 cm j 

Pp = 1.05 gm/cc “i 

Fo = 0.36 I 

Polymer dose = _J 

0.208 mg/lit PEI-18 j 


m — 





480 


Fig. 61. Model vcrificotjon Cnos 11.0 mg/lit ; dp = 0.109 jj.m) 
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Both experimental observe tiais and the model predicted 
values are depicted in Fig. 6,1 which indicate a very good 
quantitative agreement between experimental results and model 

predictions. Ihe values o£ S and R for the n/n are 1,89 and 

o 

0,99 respectively . The corresponding values for the head loss 
are 4,07 and 0,99 respectively. 

Ihe model is further verified using data from other 
experiments by Habibian ( 1971) in which one or more filtration 
variables were significantly different fjrom conditions described 
in Fig, 6,1, Observed values of the particle removal by the 
clean filters were used to calculate ri a for each run. Coeff- 
icients p / p , and p* are kept constant as they are not expected 
to vary significantly with surface properties of the media and 

particle, and concentration of particles, Ihe value of r) a 

P P 

was estimated for the best fit of the experimental results as 
it would vary vjith the surface properties of the media, suspoided 
particles, and size of the suspmded particles. 

Data from six filtration experlmoats using three suspended 
particle sizes, five influent concentrations, and two bed depths 
are presoited in Figs, 6,2 to 6,4, The model predictions are 
also shown in these figures indicating a very good quantitative 
agreemoit with the experimental results, Ihe respective 
values and the corresponding S and R values are presented in 
Table 6,2, 'Ihese values also justify the good correlation 
between experimental values and model predictions in most of 


the situations 



HEAD LOSS, cm n/Oo, percent 



80 '- 

i 

60- 




dc = 0.038 cm 

Vq = 0.136 cm /sec 1 

L =14.0 cm I 

Pp = 1.05 gm/cc 

Fo 5 0-36 I 

Polymer dose = 

0.0682 mg/nt PEI-18 ! 



TIME, min 

Fig. 6.2. Model verification (noS 9.7 mg/iitj dp s 0.109 |jum] 



HEAD LOSS, cm n/n©, percent 


60 


1“ "'T 



Experimental 
{Habibian, 1971) 
Model predicted 


7 2 




0 




120 


dr = 0.038 cm 



Fig. 6.3. Model verification (nos 37.2 mg/lit, dp=1.1 jj.m ) 



HEAD LOSS, cm n/no, percent 


“ ' * — r 


40f 


20 F 






)isiri 






'“1 '’'r 


<Jc * 0.038 cm 
Vo s 0.136 cm /sec 
L *2.0 cm 
Pp * 1.05 gm/cc 
Fo * 0.36 


▲ A 




200 r 


160!- 


120 



(1) ■, (2) A. (3) • Experimental 

" - ' Model predicted 


dp./i.m 

0.109 

1.1 

7.6 


^o, mg/ut 

50.0 

48.0 

52.0 


Polymer dose 
PEI -18 - 

mg /lit 

0.176 

0.032 

0.024 


60 80 100 
TIME, m*m 


A. A ^ A 

- - ^ m .. .w -I.- II . 

120 140 


160 


Fig. 6.4. Model vGriflcation for particle sizes in the 
different regions of filtration. ^ ^ 
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Table 6,2, Computed values of a and statistical parameters for 

ir 

the data available in literature 


Pig, 

No. 

1 

'’p 

ap 

Standard 

Estimate 

Error of 
(S) for 

Coefficient of 
Correlation ( R) for 


J {im 


1 

n/n { Head loss 
n ! 

n/n^ 

f 

; Head loss 


1 , - .. >, 

L „„ ) 



, . 

f 

6.1 

0.109 

0.041 

1.89 

4.07 

0.99 

0.99 

6.2 

. 0.109 

0.042 

2.02 

8.26 

0.96 

0.99 

6.3 

1.100 

0.140 

2.54 

2.51 

0.94 

0.99 

6.4 

0.109 

0.041 

3.25 

1.96 

0.89 

0.95 

6.4 

1.100 

0,065 

7.50 

2.33 

0.73 

0.97 

6.4 

7.600 

0.173 

0.95 

0.22 

0.97 

0.97 


n IT ^ . rt-m. 


te. - A \ r w 

r A- Air>^ -A ■ 




It is important to note that the a value increases with the 
increase in suspended particle size (Table 6.2). In other words ^ 
one may write 

m 

Up a (dp) ^ (6.11) 

Heres, is a constant. 

This of relationship of a with d is geierally expected 

ir ir 

because as the particle size decreases, the charge on the particle 
surface increases. Ihe increase in particle surface charge 
increases the thickness of the electrical double layer and hence 
the energy barrier for sucx:essful attadhamoit of the two particles 

which leads to decrease in cx value, This observed relation ^lip 
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of with suspoided particle size gives an additicsnal verifies- 
aticn of the model since it is able to cempu-ce cc values in 
accordance with the expected variation of cl with d . 

p p 

6.2 Model Verification Using the Data obtained in 

Present investigation 

The experimental investigation was carried out to suppleanent 
the data available in literature on deep bed filtration contro- 
lling the suspended particle size to check the validity of the 
proposed model under different sets of operational and systen 
variables , 

The major difference between the experimental investigaticn 
of Habibian ( 1971) and that of the present study is the use of 
alum as destabilizing agent instead of polyelectrolytes, The 
second significant difference is the initial bed porosity and 
the use of two differoit media sizes, 

Ihe data from three filtration experiments employing three 
different suspended particJLe sizes/ viz., 0,109 jim, 1,1 jim, and 
7,6 }im, are presented in Figs, 6,5 to 6,7, The model predictions 
are also shown on these figures indicating a very good quantita- 
ti\'.j agreOTient with the experimcaital results, Ihe respective 
values and the corresponding values of the statistical parameters 
are pres^ted in Table 6,3. 

It is interesting to note frc»n Figs, 6,5 to 6.7 that even 
though the initial bed porosity is same with 0,03 and 0,05 cm 
media sizes, the initial head loss, particle rotioval, and the 
head loss developnent patterns are quite differ^t. Ihis clearly 
indicates that the surface area and the individual pore sizes 
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Table 6,3, Computed values of and statistical parameters for 
the data obtained in present inve stigation 


Fig. 

No. 

d 

P 

“P 
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Estimate 
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are more important in filtraticn than the initial bed porosity , 
Furthermore, as the filtration progresses the change in surface 
area of the media and the individual pore size are mainly 
responsible in head loss development than the decrease in bed 
porosity. This is because the comjxited changes in the bed 
porosity (Eq, (4*10)) as the filter run progresses are found to 
be negligible in all the experiments. Hence there was no 
noticeable difference in model predicted filter performance with 
and without accxunting for the variation in bed porosity. Thus 
the assun^ticn of constant bed porosity is quite justified in 
the filtratiai of water and wastewater. 

To be more realistic, the model ^ould be verified for 
varying parti cdLe sizes in the influent. In order to achieve this. 
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experiments were ccncaucted with different combinations of the 
three available latex particle sizes, viz., 0.109 jim, 1.1 pm, 
and 7,6 pm, "Ihe total influent particle concentraticai was kept 
constait in all such experiments. Ihe main difficulty in these 
experiments was to estimate the effluent c^ality. Many attempts 
were made to estimate tlie concentration of different particle 
sizes in the effluent, tut none of these v;ere successful, ihe 
only method available to estimate the particle size distritutlon 
is with the use of Particle Ccaanters. This method ccaild not be 
used because of the non availability of Particle Counter, 
Therefore validity of the model could be checked cnly for the 
head loss data. The results of the fcur filtraticn experiments 
conducted with different coiribinaticns of latex particles are 
presented in Figs, 6,8 to 6,11, The model predicted head loss 
values are in good agreement with the experimental observations 
indicating the validity of the model even with varying suspended 
particle sizes in the filter influent. However, the definite 
conclusion in this regard can be drawn only when the experimental 
observations confirm the model predictions for effluent quality 
also, 

6.3 Filtration of Clay Suspension 

Filtration experiments with clay suspension using glass 
beacSsand sand as filter media were conducted with two main 
objectives 

1, to compare the filtration performaice of two media, 
viz,, glass beach and sand, and 
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2, to compare the filtration of latex suspension with 
clay suspension. 

Ihe main differences between the sand media and glass beads raedia 
were the initial bed porosity and sphericity of the particle, 

Ihe initial bed porosity of the sand bed was 0.467 while that 
of glass beadsbed was 0.429, Ihe major differerKie between the 
latex suspension and clay suspension is that the latex particles 
are exactly spherical and have a particle density of 1,05 gm/cc 
while clay particles are expected to be flaky and have a particle 
density of 2,65 gm/cc. 

Many experiments were conducted with clay suspoision using 
glass beads (0.03 am and 0,05 am) and sand (0,03 am GM) , and 
varying physical and physico-chemical parameters influencing 
filtration, Ihe results of one such typical experiment are 
presented in Fig, 6,12, Qualitatively both sand and glass beads 
behave in a similar way* Ihis indicates that the proposed model 
could as well be valid for sand media which is very commonly used 
for filtration of water and wastewater, only the mcdel parameters 
6,8, and S' would be different, similarly the <oaipa»ison of 
the filtration of clay suspension and latex particles (Figs, 

6.5 and 6,12) Shows a qualitative agreement, Ihis would mean 
that the prop>osed modei could also be applied to the filtration 
of cd. ay suspensions, 

A direct quantitative proof of the validity of the proposed 
model for sand media using day ^spensicn was necessary and 
would have be^ much more convincing to accept the mocfel. Again, 
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unfortunately the lade of particle size distributicxi data limits 
the present study to give only a qualitative validity of the 
model in such conditions. 

6,4 Variation of P arti cl e-to-P article Attadiment Coefficient 

Parti cl e-to-particle attachment coefficient (oCp) would 

mainly depend upon the surface properties/ particularly the 

surface charge of the particles. The particle surface charge 

may be expressed in terms of 2eta potential of the particles. 

In other wordS/ a can be related to zeta potential. In coagu- 

lation and filtration the zeta potential is controlled by the 

addition of destabilizing agents such as bi- and multi- valent 

cations, poly electrolytes etc. It would be very much useful if 

a relationship between a and the zeta potential can be developed. 

P 

However, non existence of sufficient data dxxit filtration 

effici^cy controlling both particle size and zeta potential 

limits the development of such a quantitative relationship. 

Presently an attempt has been made to ei:ipirically correlate 

a with poly electrolyte dose which in turn controls the zeta 
P 

potential, Ihe data of Habibian ( 1971) have beoi used for this 
puriose. 

The computed a values aire plotted with polyelectrx)lyte 
P 

dosages for two different particle sizes in Fig, 6*13. Noting 
the nature of the variation of with polyelectrolyte dosage 
an empirical relaticn is pixjposed relating these two variables. 
Ihe equation is of the following form 

. . 2 

Ln et = - (A. X + A- X + A_) (6,12) 

p 1 2 3 



Ln {a. 
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Hera, x = Ln (Polymer Dosage, mg/lit), and and are 

empirical constants. 

Obviously, one \(C>uld expect A^, and to depend upon the 
particle concentration, suspended particle diameter and surface 
charge of the particle. 

Fig, 6,13 clearly indicates that a is more soisitive to 

P 

the polymer dosage for suspended particle size of 0,109 |im 
compared to that of the 1,1 jim. Hence it is more important to 
control the polymer dosage when filtering a suspension of lower 
particle size, i.e., sutmicron particles compared to the filtra- 
tion of suspension having higher particle size. Ihus the 
knowledge of particle size can help even in controlling the 
dosages of destabilizing agents. 
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7, SIMULATION 

Sirnulation assumes that the systori can be described in 
terms acceptable to a compjting systen. In this regard, a key 
concept is that of a "system state description". If a astern 
can bo characterized a set of variables, with each cxmbin- 
ation of variable values representing a unique state or condition 
of the syston, then manipulatjon of the variable values simulates 
movanent of the systan from state to state. . Precisely the 
simulation is the representation of the dynamic behaviour of 
the system by moving it from state to state in accordance with 
Well defined operating rules. 

With the aforementioned conc^t of simulation, the deep 
bed filtration can be simulated by experimenting on a verified 
proposed model for different ccambinations of varims physical 
and physicxt-chemical parameters influeicing water and wastewater 
filtration. In general, the simulation procedure involves the 
prediction of the dynamic behaviour of 

1, the system performance, and 

2, the system state description. 

In the present context, the ^stam performance refers to 
the filter performance whicdi includes filtrate quality and head 
loss* He systCTi state descripticn refers to the state or 
condition of tha filter bed which in turn gives the clogging 
picture of the filter l:^d. The simulation of filter performance 
cffid the state of the filter bed are descrribed in two separate 


secticxis 
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7^1 Simulation of the Filter Performance 

The variables influencing filtration may be classified in 
two groups / viz* physical parameters and physico-chemical para- 
meters* The physical parameters include suspended particle 
concentration/ density and size/ media size/ porosity/ flow ratO/ 
and filter bed depth* The phy si co— chemical parameters include 
particle— to— media grain attachment and particle— to— particle 
attachment* The simulation of filter performance for all these 
parameters is described as follows 

A, Su pended Particle Conc^tration: Hie predictions of 
effluent quality and head loss fer four different influent 
suspended particle concoitration are presented in Fig. 7,1, At 
the beginning of the filtration, removal efficiency and head loss 
are independent of particle cxjnc^tration as it is not expected 
to play any role initially, and hence properties of the filter 
bed and suspended particles will govern the filtration and not 
the concentration of particles in influent. However, the 
dynamics of filtration completely changes with the particle 
concentration. As particle cmceitration increases removal 
efficiency increases because of the increasing number of retained 
particles acting as collectors. After some time the filtration 
efficiency attains a saturation value because of the surface 
coverage of the media grains and retained particles acting as 
collectors, A dynamic equilibrium may be achieved as the new 
contact sites are not effectively increased and is influehoed 
the particle concentration. As expected, heai loss development 
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removal at the onset of filtration. This is in accxrrdance with 
the observations of Yao Qt.al, (1971), Further/ CTie micron 
particles are the most difficult to ronove as the filter run 
progiresses. This is because the rate of increase of retained 
particles acting as collectors is very slow when only 1.0 micron 
particles are present. Submicron particles give better effi- 
ciency as th^ can be effectively transported by Brownian 
diffusion. These particles produce very hi^ head loss because 
of their more surface ’area per unit weight of the material removed 
as compared to bigger particle sizes. The bigger particles 
(5.0 |i.m) show more removal as they can be very easily transported 
by interception and gravity settling. Further, because these 
particles have low surface area per unit weight of the material 
r^noved, the head loss developmoit is not significant. This is 
in accordance with the experimental observations of Habibian 
(1971), and that of the prescaat investigation. 

D, Varying Suspended Particle Size in Influent: In order 
to simulate the filter performance in more realistic cases, the 
influait is considered to have three different particle sizes 
nanely 0,109 ^m, 1,1 ^m and 7.6 jim. The mass cxjncmtration of 
all these particles is tak^ as same. It is further assumed 
that the surface properties of all these particdLes are such that 
they have equal opportunities for aiccessful contacts. The 
results of simulated filter performance are presented in Fig, 7,4, 
The removal of 1,1 p.m particles Continues to pose a problo^i. 

The effluQit 'mDuld contain more number of 1.1 ^im particles 



HEAD LOSS^ cm n/no, percent 


100 


40!- 


20^ 


97 





t. ,. a Jtir 


^ w u 


I r r ikf 


1 

4 




jn ri 


0 

50 

40 

30 


20^ 


10 p 

i 




1 

2 

3 

4 



Jn»*. .* . ■ ^ * 1 ■<» r /■> 1 r n » •■ Mira ■ 


Hi /Ho, 

Oo, (dp^ = 0.109 fi,m) = 6.65 mg/lit 

-i- » *.>3 

1 

? 


nojCdpj = 1.1 |j/m) = 6.70 mg/lit 



nojCdpjS 7.6 ji/m) = 6.65 mg/lit 

1 

£ 

n /no j 

Oo - Oo^ no 2 ^03 

J 

{ 

* 

n s n| +n 2 + n 3 

0 

1 

\ 

i 

i 


’ no = 20.0 mg/ lit 

i 

1 

fc 


dc = 0.038 cm 

J 


Vo = 0.136 cm /sec 

1 


L = 15.0 cm 



Pp = 2.65 gm/cc 

' 


Fo:= 0.36 

5 


^0 


... ,J,./ 

40 


L„ 

280 




320 


I, .. . . . 1 . . 

ii f s 120^ 160 200 240 

Illllii 

Fig, 7-4. influence of varying particle sizes in suspension 
oft^' liltir^ = performance - ■ 



98 


comparod to 0.1 [im particles. Hie 7,6 fm particles would be 
removed almost completely . in general it may be said that the 
filtered effluent would always have major fraction of particles 
around 1,0 p,m size. The head loss development is a direct 
consequence of the extent and the type of particles ronoved. 

E. Media Size: Computations illustrating the effects of 
media size on filtration are presented in Fig, 7.5, As expected, 
both removal efficiency and head loss increase with the decrease 
in media size initially as well as with time, 'Ihcugh, media 
size affects the filter performance, it can be easily measured 
and readily controlled, and haice is not of major importance, 

F. Clean Bed Porosity: The simulated filter performance 
for four differ^t clean bed porosities is presented in Fig, 7,6. 
It is interesting to note that tdie head loss baild-up increases 
drastically with decreased clean bed porosity whereas the ranoval 
efficiency is not affected to that extent, This is essentially 

■4 

because of the more uniform distribution of particles with depth 
(Refer Part C, Section 7,2). 

G. Flow Rate: Flow rate is an important operational 
variable to be contrxDlled in filtration process. As is clear 
from Fig, 7,7, flow rate affects both the initial filtrate 
quality and head loss. However, as the filtration progresses 
the roaoval efficiency becones independent of flow rate , The 
head loss development, of course, is a direct consequence of the 
cumulative depositica:! in the bed and the flow velocity. 
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Pi. Filter Depth: The effects of filter d^th and filtration 
time are illustrated by the simulation results presoited in Fig. 
7.8, Removal is distributed throughout the filter bed at the 
onset of filtration, hut becomes localized in the upper regions 
of the bed as filtration proceeds. Head loss follows a similar 
distribution. Retained particles accumulate in the upper regions 
of the bed and lead to the retention of more particles in that 
zone, A good filter is one which would give uniform rsnoval and 
head loss with depth. This can be achieved by properly control- 
ling the system and operational variables in filtration. For 
example, the uniform ronoval with depth can be achieved by provi- 
ding larger sized media with lower specific gravity at the top 
and smaller media with higher specific gravity at the bottom 
(mixed media) or by adopting upflow filtration. Ihese observa- 
tions cJL early justify the superiority of multi media filters over 
single media filters. The rate at which the retained particles 
would increase is dependent on the m«3ia size. If the media size 
is large at the top it will give lower removal at the txjp and 
hence the particles are forced to poietrate deep whore the media 
size is low. At lower d^ths, thou^ the media size is low, 
because of the low particle concentration at that depth, the rate 
of accumul ation of retained particles is also low, Thus by 
adjusting the media size*/ rsmoval can be distributed throu^iout 
the d^th. Further the rartK>val efficiency will v^uy with other 
system and operational variables and can be predicted from the 
proposed model which in tu^ can help in manipulating these 
Variables dep^ding upcsi ^steta input variables 
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a, and a^) and the output requirements. 

I. ParticlQ-to-Media Grain Attachment: Ihe effect of 
particle- to-media grain attachment on filter performance is 
depicted in Fig. 7.9. This is adiieved by varying the coefficient 
ct in a simulation model which is a reflection of the particle- 
to— media grain attachment. Initially the r^oval increases with 
increase in a value. However/ as the filtration progresses, the 
removal efficiency might increase or decrease depending upon the 
^p Hiis is because as the filtration progresses the media 

grain gets covered due to deposition of particles and the proper- 
ties of the deposited particles control the attachment step. 
Depending upon the relative magnitude of the decrease in value 

due to surface coverage and the increase in ti value due to the 

r 

retained particles acting as collectors the overall removal 
efficiencY may increase or decrease. The particie-to-media grain 
attachment (a) is an important parameter and can be controlled 
by controlling the media and particle properties hy giving proper 
pretreatment either to media or susp^ded particles or both. 

This is normally achieved by giving precoat to the media and 
adding poly electrolytes to the suspensions or by giving different 
pretreatmoits to the filter influait. Ihe major facrtors which 
are expected to have influence aa a are the zeta potentials of 
media and suspended particles, 

J, Particl#-to-^artic^ P article- to-partica.e 

attachment plays a dominant irole in filtraticai as the filter 

run progresses, Th4f|, effect: = is studied by varying the coefficient 

*ct * in a sijmilaticn moclel. ; results are presented xn 
P ‘ 
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Fig. 7.10. If oip is zero, i.e., there is no particle^to-particle 
attachment, then the filter efficiency wjuld continuously 
decrease because no retained particles would act as collectors 
and at the same time these retained particles would cover the 

3-^®^ ®f the media thereby decrease the particle— to— media 
grain attachment and hence the overall removal efficiency. As 
oCp increases, the removal efficiency as well as the head loss 
increase significantly with time. On the other hand, the particle— 
to— particle attachment does not have any effect on the removal 
efficiency of clean filters. 

7,2 Simuliation of the state of Filter Bed 

The retention or deposition of the suspended particles as 
the filtration progresses results in the change of the condition 
or state of the filter bed These changes are indirectly reflected 
by the output variations characterised ty the filter performance. 

In order to optimize the filter response for the different sets 
of input conditions it is necessary to choose the appropriate 
values of the operational and system variables. In this respect 
the dynamic changes of the state of the filter bed can play an 
important role. The filter state may be characterized by the 
clogging zone. This corresponds to that portion of the filter bed 
which contributes for the head loss tuild-up. The depth of this 
zone changes as the filtration .progresses, and at any time ( t) 
after the start of filtration the depth of clogging zone is 
governed ly the depth of floe p^etration or the d^th of ciLogging 
front. The concept of clogging f rent, was first introduced by 
Stanley (1955), and latter on ^ployed by several investigators 
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(Crapps, 1964; Huason, 1964; Adin and Rebhun, 1974). However, the 
theoretical formulation to obtain the clogging front is not yet 
available. 

The clogging front is defined by the points where the head 
loss gradient begins to deviate significantly from the value 
for a clean filter bed. Mathematically the clogging front may 
be defined by tbe locus of the point where the rate of change of 
the head loss gradient beccmes zero^ and may be ej^ressed the 
following equation, 

^(dhf/dL) = O (7.1) 

or 

,2 

*^'“2 ^ ^f ^ (7.2) 

dL 

Thus using Eqs, (4,7) and (7,2) the dynamic behaviour of 
the clogging front which would indirectly define the filter state 
can be predicted. The simulation of the dynamics of the clogging 
front for various physical and jhysico-chemical paran^ters is 
described as follows. 

A. Physical Properties of the influent Suspension: influence 
of suspended particle concentration, density and size on the 
variation of the depth of clogging front and the head loss at the 
depth of clogging front are depicted in Figs, 7,11 to 7,13 respec- 
tively, ihe depth of clogging front increases with the increase 
in suspended particle concentration and decrease in particle 
daisity (Figs, 7,11 ^ ‘7^12)*: HcfWever, the head lo^ at the 
depth of clpggii^ frrm^ a^fpairas w^aitered either with the variation 
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in suspended particle concentration or density. Here again the 
su^ended particle size show a typical behaviour, one micron 
particles take tlie longest time to form a noticeable clogging 
zone, Sutaaicron particles form a deeper clogging zone from the 
beginning of the filtration process, Ihe larger particles also 
form the clopging zone from the beginning but the depth is very 
very small, 

B, Media Sizes Media size plays an important role in the 
penetration of the floes. As the media size increases the depth 
of floe penetration increases. Further the head loss build— up 
at the d^th of clogging front is lower for hi^er media sizes 
(Fig, 7,14), The size of the media can be controlled very easily 

and hence depending upon the input conditions (n / p , d , a, and 

^ 

ffp) the appropriate media size can be choosen to achieve the 
required filitrate quality and more uniform distribution of the 
particle by predicting the filter performance and the nature of 
clogging front using the proposed model, 

C. Clean Bed Porosity s The variation of the depth of 
clogging front and head loss at the d^th of clogging front for 
four different initial bed porosities are plotted in Fig, 7,15. 
The depth of floe penetration increases with the decrease in 
initial bed porosity, ihe head Ib^s at the depth of clogging 
front is lQ,wer for lower depths of ciLogging front with increase 
in bed porosity while this t^d is reversed for hi^er depths of 

clogging; f.pht.:.iF’ii'i'^'|k|§fc^iiilii|iij^ ' ’ - 
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D. Flow Rate: As expected the depth of clogging front 
increases with the increase in flpw rate (Fig, 7,16). However, 
the head loss built at the depth of dogging front follows a very 
typical trend. For lower depths of clogging front the head loss 
is higher with higher flow rates whereas at higher depths of 
clogging front the head loss is higher with lower flow rates. 

This clearly indicates that with higher flow rates the deposition 
is more uniform, 

E. Physico-chonical Parameters: The physico-chemical para- 
meters are represented by partide-to-media grain attachment aid 
particle— to-part id e attachment. The variation in the depth of 
clogging front and its relationship with head loss for physico- 
chemical parameters may be studied by varying the coeffici^ts a 
and (Xp in the proposed model. The simulation results are presented 
in Figs. 7,17 and 7,18. The partide-to-media grain attachment 
does not have much significance with respect to the depth of floe 
penetration. The control of particle-to-partide attachment to 
seme extent can distribute the removal over the d^th of the 
filter bed. The head loss at the d<^th of dogging front is 

of course a direct consequence of the amount of deposited material 
and in general has no relation with either a or 
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Fig. 7.17. Effect of particle -to-filtcr grain attachment on 

temporal variation of the depth of clogging front 
and its. relationship with head loss. 




Fig. 7.18. Effect of particle -to -particle attachment on tempoml 
variation of the depth of clogging front and its 
relationship with head loss. 
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8. SUMMARY 

In the pjTesent study a sirrajlation ptocedure is developed 
for the dynamic predicticxi of filter performance and filter state. 
The method is based on ^nthesizing available q[uantitative results 
relating to filtration and to porous media flows, within an over- 
all fraraevtork which views the process to consist of two principle 
st^s. namely the transport of the susp^ded particles to filter 
media grahis or already retained particles and the attachment of 
these particles to the filter media grain or retained particles 
acting as collectors. Evaluation of the results through compa- 
rison with the data available in literature and the data obtained 
bhe ej^erimental investigation c£ the presoit study indicates, 
to the extent the nature of these types of data permits, that the 
method is very effective and even capable of predicting on satis- 
factorily quantitative basis some intricate details of observed 
filter behaviour. 

In the simulation procedare, the filter performance is 
characterized by the filtrate quality and the head loss develop- 
ment, The filter state is defined by the depth of the clogging 
front which decides the clogging zone in the filter bed. Ihe 
results of the simulation studies indicate that the proposed model 
is able to predict the filter performance in accordance with the 
observations made in laboratory and field. Based on the results 
of the simulation studies following significant obseirvations are 
made which may lead to the process modification and/or development 
of suitable design criteria for the deep bed filters. 
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Ttie major important input parameters which affect the 
filter performance are the size of the suspended particles and 
their distribution. As the suspended particle size decreases the 
head loss development increases for the same OTount of r^oval , 
The filtration efficiency thou^ influenced by particle size^ a 
goieralized statonmt regarding the variation of filter efficiency 
with particle size is not possible because it is affected ly both 
transport and attachment step. The transport efficiency is 
minimum for depended particles around one micron while the 
attachment efficiency is expected to increase with the increase 
in particle size. Ihe overall removal efficioicy depends upon 
the resultant effect of both transport and attachment step. 

The another important input parameter is the density of the 
suspended particle which controlles the raimber of particles in 
the suspension if the filter influent concoitration is expressed 
in terms of mass per unit volume. To be very precise the filter 
load should be expressed in terms of the number of particles per 
unit volume with particle size distribution rather than expressing 
the same in mass per unit volume or turbidity units. 

There are several other system and operating variables which 
influence the filter performance whose effect can vesy well be 
studied from the proposed model and can be manipulated in the 
actual practice to get the optimum output conditions with respect 
to filtrate quality and head loss development. The depth of 
clogging front may be used as a criteria for the design of filter 
depth. 
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9 , RECOMMENDATIONS FOR iUTURE RESEARCH 


On the basis of the results of the csirrent investigation 
the following reccanmendations are made for the future research as 
a logical continuation of the research presen ta3 in this thesis. 

1. The, pressed model should be verified for varying suspended 

particle in the filter influoit. 

2, Ihe simulaiticai model should be modified to incorporate the 
reentrainm<a3t of the deposited particles. 



possible chcttiges in the size of the suspoided particle 
within the pores of the filter be3 beeause of coagulaticn 
should be predicted and appropriately incorjx) rated in the 
simulaticai model, 

4, The model should be extended to include the effect of 
irregular shapes of the suspended particles and the media 
grains. 

5, A better understanding of the influence of physico-chonical 
parameters on filtration is necessary. This may be done 
utilizing the knowledge of coagulation and stating the 
similarities between coagulation and filtration, A model 
should be formulated relating a and to surface charac- 
teristics ( zeta potentials) of the media grains and 
suspended particles. This warrants detailed experimental 
investigation on coagulation and filtration with the 
proper control of the surface properties/ in particular 
surface potentials of the media and suspended particles. 
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Hae aforementioned in^^'ovanentS/ however, cannot be made 
with further experim^ts ^rfonned under well defined conditions. 




The importance of d^ta gathered from a systematic programme of 
critical self-evideit. . 
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APPENDIX A 

SOLUTION TECHNIQUE 

Equations (4,2), (4,3), and (4.4) constitute the mathe- 
matical model for predicting the filter effluent quality . 
Ihese equation? can be rewritten as 

N dp/dl + N r,c<p(dj/d^) = 


+ (1 - N 


(A.l) 


IN . , . IN 
3t o 3L 


= [C^ + N(C 3 - C 3 ) - N C^] n 


(A. 2) 


-~ + V Iy +'-4 “ V ^ ^ 

3 t o?L 2 d o r 


(A,3) 


Ihe analytical solution of the above equations is not 

possible. Therefore, simplifying techniques are used to obtain 

a reasonably accurate solution. Ihe algorithm to solve these 

equation proceeds as follov7s. 

Considering n and N as step functions of time and length 

th 

instead of continuous ones, Eq. (A.l) for the i time interval 
"tin 

and j length element may be written as 

d 2 

\ = na (a; - Nj ap)/d^ + p^ Nj (3^) 

j ,1 


d 2 


J. (1 - P^) Up o:p (3^) 


(A. 4) 
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with the boundary condition 


Ni^l to 


^ - = 0 resulting in n to n = r]a 

' * : 1,1 NL,i 


Here, NL is the total number of length elemoits. 

Assuming h to remain constant for a very small time 
interval /it, Eq, (A, 3) for the i time interval reduces to 


V v n. TL 

o 9L 2 d o 1 r, 
c X 


= O 


(A.5) 


Ihe above equation for the j length element may be written as 


3 ( 1 - F) 

*"^,1 = - 2 


(A.5) 


|p;ere, &L is the thickness of the j element. Further one may 
write 


n, ... = Ti - .+ & n J . 

.3+l,i J/i J/i 


(A. 6) 


with the boundary condition 


"l.l "l.NT = "o 

Here, NT is the total number of time intervals. 

Similarly, Eq. (A. 2) for the i^ time interval and 
length step may be written as 


N, . , = N. + AN. . 

■j,i+l 3/i 3/i 


(A. 8) 


with the initial boundary condition 
^'^1,3 ^1,NL ^ 

A flow sheet illustrating the above procedure is given in Fig. A.l. 
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No^ Noi 




^ Specify the time interval At I 
: and length step AL I 

N = n 

Ujli 


n — ni 


.«.* I 

\ N -Nj,i J 

-*? z ^rj.i . \ 

i [Eq. (A.4)J I 
j TL r 0.0 I 

nj.i : nj,i ♦Anj.i ^ 

[Eqs. (A. 5) and {A.6)j 
Nj.ui =Nj.i.ANi 
[Eqs. {A.7)and(A.8)j 


I j = j 

1 TL =TL * AL 


-{is TL : Filter depths 

„ JYgs ^ 

\ Time = Time + At | 

; i = i + 1 i 

- • 

^ Is Time > Filter run . 


lYcs 



t A.I. Proposed (ALGORITHM) logic diagram for predicting 
filtrate quality . 
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The procedure used was proposed by Marquardt ( 1963) as 
an extension of the Gauss Newton method to allow for convergenGe 
with relatively poor starting guesses for the unknwon coeffi- 
cients. In this method, the Gauss Newton normal equations are 


modified by adding a factor X/ 

^ i . I] iA = A^(y_”y*) 

Here, 



(B.3) 




13.3 


I is the identity matiriX/ A is the transpose of A matrix, 

and the asterisk designates quantities evaluated at the initial 
trial values. 

Thus X is added to eacii term of the main diagonal of the 
A _A matrix, Ihe rules for calculating ^ are discussed in the 
original article ty Marquardt ( 1963) . All other compitaticaial 
details parallel the Gauss Newton procedure. 

The ccntputer progrartme utilizing aforementioned algorithm 
is available in "Optimization Techniques with Fortran" (Kuester 
and Mize, 1973) , The progranrae cxinsists of a main prograimie, 
a general subroutine BSCajVE, a general functim subprograrrane 
Sii?COSS, and tv7o user supplied subroutine^, FUNG and DERIV. 

All input and output is through the main programme, ' 
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appendix c 


STATISTICAL PAR^^IETERS 




T^e important statistical parameters which are 

generally used to check the accuracy of idle models are 

( 1 ) CoefficiQit of Correlation 
< 2 .) Standard Error of Estimate 

Hiough the coefficient of correlation has limited 
^licabiliiy in the nonlinear case (Tukey, 1977) it was used 

only for comparative purposes along with means and Standard 
Error of Estimate, 

( 1) Expression for coefficient of correlation (R); 


n 


R 


.1 ^i^i " n X y 
= 1=1 ,, . 1 


s s 

X y 


n - 1 


(c.l) 


H^re'i 


; 

n 


2 


V 


I “ X) 

i=l 1 
' n‘ 1 


n 


' .1 

i=l 1 


y) 


n » 1 


.th 


= i experimental value, 

. th ^ ^ 

^i ~ 1 model predicted value , 


X 


= average of the observed values. 
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y = average of the model predicted values, and 
n = total number of observatiOTis , 

( 2) Expression for standard er^gr of estimate ( S) : 


Here, 

. th . 

= 1 experimental value, 

= i model predicted value, 

n =: total number of observations, and 
K = number of unknovns. 


I' 


(C.2) 



